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I . INTRODUCTION 

V 

This  report  contains  selected  documentation  of  the  MI COM  real-time 
STINGER  simulation.  Sections  II  and  III  describe  the  MAJN  program  which 
serves  as  a driver  for  pre-  and  post -real -time  computation,  respectively. 

Other  codes  which  support  these  computations  are  presented  in  Sections  VI, 

VII,  VIII,  IX.  The  real-time  computations  are  performed  by  Subroutines  FLIGHT 
and  REALT  which  are  documented  in  Section  IV.  The  interconnection  of  various 
computation  elements  is  described  in  Section  V which  is  titled  Real-Time  I/O. 

In  Appendix  A the  aforementioned  code  is  presented.  Then,  in 
Appendix  B a restructured  version  is  given.  The  restructured  documentation 
can  be  of  considerable  value.  However,  this  value  will  only  be  realized  from 
future  development . 


1.0 


II 


PROGRAM  MAIN  (PRE-REAL-TIME) 


±3&JE£. 


DYNAMICS.  INC. 


II.  PROGRAM  MAIN  (Pre -Real -Time) 


f 


Program  MAIN  provides  both  pre-  and  post -real -time  functions. 
However,  for  purpose  of  documentation,  program  MAIN  will  be  treated  as 
two  separate  pieces  of  code  dealing  with  pre-  and  post -real -time 
functions.  In  this  section,  pre-real-time  computations  performed  by 
program  MAIN  are  documented.  Portions  of  program  MAIN  not  included  in 
this  section  are  discussed  in  the  next  section. 


Functionally,  program  MAIN  serves  to  connect  on  a macroscopic 
level  several  physically  separateable  modules  of  code.  In  some  cases 
however,  computer  implementation  does  not  lead  to  a strictly  modular 
code.  Therefore,  some  portions  of  program  MAIN  involve  the  micro -codes 
which  are,  in  concept,  portions  of  other  major  modules. 


Non-real  time  data  input  to  the  STINGER  digital  computer  program 
is  organized  into  tiiree  distinct  levels.  These  three  levels  of  data  input 
are  related  to  frequency  of  data  change  during  normal  simulation  activities. 
Specifically,  data  can  be  input  to  the  program  in  the  following  ways. 


• Data  card  or  DDS  commands 

• Block  Data 

• Inline  computer  code 


Data  input  via  block  data  or  inline  computer  code  requires  careful 
attention  to  the  particular  code  affected.  Pre-programmed  diagnostics  and 
informative  printout  will  not  generally  reveal  errors  in  such  modifications . 
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More  information  concerning  pre -real  time  data  processing  is 
contained  in  the  Sections  titled  SUBROUTINE  KSCALE,  INIT,  and  GUID. 


FIGURE  2.1  DATA  CARD  OR  DDS  DATA  INPUT 


2.1 


TARGET  ACCELERATION 


= Wi+1) 

Y^(i)  = Y^Ci+l) 

i^(i)  = ij^Ci+1) 


Xpg(i) , i = 1,  27 

Y^Ci),  i = 1,  27 

Y^Ci),  i - 1,  27 


NOTE:  AT,  the  time  between  samples  in  the 
target  velocity  tables,  is  assumed 
to  be  one  second. 


TEST  FOR  STATIONARY  TARGETS 

ww  = Z7!  *7  *7 

XfgCO)  + YfgCO)  4 Z^CO) 


WW  > 338  FT/SEC  -*■  moving  target 
WW  < 338  FT/SEC  -*■  stationary  target 

TARGET  ALTITUDE 


h T = h + Ah 
asL 

where  Ah  is  altitude  of  launch  point  above  sea  level 


L 


2.5 


A 


Target  velocity-time  maneuver  tables  are  shifted  in  time  by  the 
parameter  TAU.  Setting  TAU  to  a positive  real  value  causes  the  target 
trajectory  time  = TAU  and  real  time  = 0 to  coincide.  Therefore  the  portion 
of  the  target  maneuver  tables  between  time  = 0 and  time  = TAU  is  disregarded 
The  end  of  the  target  maneuver  tables  is  filled  with  a constant  velocity 
maneuver,  i.e.,  the  portion  of  the  trajectory  between  time  = 27-TAU  and 
time  = 27  seconds.  For  TAU  = 0 there  is  no  shift  in  the  target  velocity  - 
time  maneuver  tables. 


EXAMPLE:  TAU  =2.5 


TM  (1) 

XDEM  (J , ITT) 

TMA(I) 

XDM(I) 

0 

XDEM (1 , ITT) 

-.5 

XDEM  (3 

1 

XDEM (2, ITT) 

.5 

XDEM  (4 

2 

XDEM (3, ITT) 

1.5 

XDEM (5 

23 

XDEM (23, ITT) 

22.5 

XDEM  (2t 

24 

XDEM  (24, ITT) 

23.5 

XDEM(2< 

25 

XDEM (25, ITT) 

25.0 

XDEM (2( 

26 

XDEM (26, ITT) 

26.0 

XDEM (2( 

■*-  t = 0 


LAUNCH  AZIMUTH  AND  DIVE  ANGLE 


Calculation  of  launch  angle  and  dive  angle  are  based  upon  the 
following  definitions  and  conventions: 


V - Projection  of  target  initial  velocity  vector  on  to  the 


plane 


W - Projection  of  initial  Lqs  vector  in  the  X^  - 
6 - Angle  from  positive  X„  - axis  to  V ' 


plane 

E (positive  clockwise) 

Zo-  Crossing  angle  = angle  from  V to  W (positive  clockwise) 

i|>L~  The  angle  from  positive  X - axis  to  the  vector  from  the  origin  to 

. ^ T7\ 


target  (note  this  is  just  -W) 


From  the  above  definitions  Zo  + B is  the  angle  from  positive  Xe  - axis 
to  vector  W.  Then,  since  -W  is  in  opposite  direction  from  W,  the  angle  differs 
from  Z0  + 8 by  180  degrees.  Therefore 


<|>  = Z + 8 ± 180° 

L 0 

Here  is  chosen  base  on  whichever  sign  yields  a positive  value. 

Based  upon  the  value  of  X^g)  and  Y^g)  the  value  of  t^O)  is 
determined  via  four  special  computational  rules.  Note  that  the  case 
XTL  (°)  ¥ 0 an*3  YrpL  (°)  = 0 is  excluded  from  consideration  and  should 
not  appear  in  candidate  trajectories.  The  four  computational  rules 
as  they  are  categorized  in  the  computer  code  are: 


CASE  1 


1 

0 

S >o 

8 = tan-1  (X  (0)/Y  (0) ) 

TE  TE 


\J;  = Z + 6 - 180°,  i|j  (0)  is  determined  by  whichever 

L 0 L 

Z^  + 8 + 180°  quantity  is  positive 
SEO  = (see  computation  in  Subroutine  INIT) 


CASE  2 


Xm  (0) 

Y (0) 

TE 

TE 

0 

>0 

3 = 90° 

tp  (0)  = E + 8 - 180°  , ip  (0)  is  determined  by  whichever 
L 0 L 

Eg  + 3 180°  quantity  is  positive 

SEO  = (see  computation  in  Subroutine  INIT) 


CASE  3 


X (0) 

TE 

Y (0) 
TE 

>0 

>0 

>0 

<0 

<0 

>0 

<0 

<0  }Cc 

6 = 270 


ip  (0)  = E + 3 - 180°  ip  (0)  is  determined  by  whichever 


EQ  + 3 + 180 


quantity  is  positive 
SEO  = (see  computation  in  Subroutine  INIT) 


CASE  4 


X (0) 
TE 


0 

>0 

>0 


Y (0) 
TE 


and  Z (0) 
TE 


0 

® ^Computer  implementation  does  not  take 
®^care  of  these  two  cases 


♦_  (0)  = 0 
Li 


For  stationary  target  crossing  angle 
is  not  applicable,  program  sets  KKK  = 1 


CASE  4b 


CASE  4c 


X (0) 
TE 


>0 

<0 

vo) 

= 0 

SEO  = 

90° 

YrPT,(°} 

TE 


0}Computer  implementation  does  not 
0}take  care  of  these  two  cases 


For  vertical  target  trajectory  crossing 
angle  is  not  applicable,  program  sets  KKK 


X (0) 

Y (0)  and  Z < 0 

TE 

TE  TE 

0 

0 

>0 

0} Computer  implementation  does  not 

<0 

0}take  care  of  these  two  cases 

V°) 

= 0 

For  vertical  target  trajectory 
crossing  angle  is  not  applicable 

SEO  = 

-90° 

program  sets  KKK  = 1 

TRANSFORMATION  OF  TARGET  VELOCITY  COMPONENTS  FROM  EARTH  FIXED  TO  LAUNCH  COORDINATES 


X_  (0> 
TL 


(0) 

TL 


2 (0) 

TL 


cos  0 cos  lj;  cos  0 sin  \J/ 

L L L L 


-sin  ifi 


cos  i li 


sin  0 cos  <1j  sin  0 sin  \l> 

L L L L 


where  0 = 0 (0) 

L L» 


*L  = V0) 


— 

— — 

-sin  0 

L 

X (°) 

TE 

0 

Y__  (0) 

TE 

cos  0 

L 

2^  (0) 

TE 

.»  - 

This  transformation  is  performed  in  Subroutine  GUID 


RANGE  TABLE  GENERATION 


1 


The  range-table  is  generated  for  use  as  a data  collection  tool 
during  real  time  simulation.  The  range  table  is  a set  of  relative  missile 
minus  target  positions  along  the  XL  axis  at  which  data  collection  should 
be  performed.  Generation  of  the  range  table  follows  from  specification  of 

XTERM,  the  last  data  point  to  be  collected 
DX,  the  distance  between  data  points 
NPTS , number  of  data  points  desired 


Provided  analog  scaling  permits  the  collection  of  NPTS  points , 
this  is  done.  Otherwise,  the  program  automatically  reduces  the  value  of 
NPTS  until  analog  scaling  requirements  are  satisfied.  The  new  value  of 
NPTS  is  stored  for  later  usage  as  NADJ. 

The  rule  for  reducing  NPTS  to  a final  value  of  NADJ  is  given 

below 


NADJ  = NPTS,  if  DIF  > 0 

NPTS  - truncated  value  of  (DIF/DX)-1,  if  DIF  <_  0 

Here  DIF  is  defined  as  the  difference  between  the  data  collection  distance 
commanded  and  the  largest  data  collection  distance  possible  with  the  analog 
scaling  given.  DIF  can  be  computed  as  follows 

DIF  = {XTERM  + (NPTS-l)DX}  - {(Analog  Scale  Factor)  (Largest  Voltage)} 

1 + Y 


2.13 
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TABLE  2.1  DESCRIPTION  OF  MAIN  PROGRAM  SYMBOLS  (Pre -Real -Time) 


PROGRAM  SYMBOL 


VARIABLE 


DESCRIPTION 


Print  flag,  IS1  ? 1 indicates  print  is 
desired. 


NSAM* 


ICON* 


CLMS(I)* 

CLS1(I)* 

CLS2(I)* 

Ql,  Q2,  ...Q6* 
Wl,  W2,  ...W6* 
NRUNS* 

NPTS* 


XTERM* 


ZTEO(I)* 


ASC(I)1 


Quality  control  test  case  sample  size 
2 <_  NSAM  <11 

Automatic  or  manual  launch  set  selection 
(ICON  = 1 implies  automatic,  ICON  = 2 
implies  manual) 

Number  of  first  launch  set  desired  in 
automatic  mode,  1 < INL  < 60 


Used  in  computation  of  quality  control 
limits.  For  further  information  see 
SUBROUTINE  QCLIM 


Number  of  runs  (NRUNS  <_  20) 

Number  of  data  points  to  be  collected 
in  real  time 

Distance  between  real-time  data 
collection  points  (FT) 

Distance  between  target  and  missile 
at  last  data  collection  point  (FT) 

NL  + 1 is  number  of  launch  sets  to  be 
read  from  data  cards,  if  NL  >_  61 

Launch  set  target  altitude  in  earth 
fixed  coordinates  (FT)  (Note:  printout 
shows  - ZTEO(I)  ) 

Launch  set  target  crossing  angle  table^ 
measured  positive  clockwise  from  V to  W 
(DEG) 


*Data  input  from  card  reader 
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PROGRAM  SYMBOL 

VARIABLE 

DESCRIPTION 

RR(I)  * 

% 

Launch  set  target  initial  range  table 
(FT)  (distance  measured  along  the  line 
of  sight) 

III* 

Flag  for  automatic  selection  of  target 
trajectory  III  = 1,  for  manual  selection 
set  III  / 1. 

IRSS* 

Flag,  IRSS  simulation  requires  IRSS  = 1, 
for  MI COM  Hybrid  set  IRSS  / 1. 

KA* 

Index  of  target  launch  set.  This  is 
used  in  manual  mode,  i.e.,  Ill  = 1 

ITT* 

Trajectory  number,  ITT  = 1 selects  a 
stationary  target.  In  general 

1 1 ITr  i 45 

N OOO* 

TAU* 

T 

A parameter  which  selects  a portion  of 
the  target  velocity-time  manuever  table 
for  use  in  the  simulation.  This  para- 
meter shifts  time  forward  TAU  seconds 
and  disregards  all  manuever  information 
between  time  = 0 and  time  = TAU.  Note: 

0 < TAU  < 27.  With  TAU  = 0 the  manuever 
tables  are  unchanged.  (SEC) 

ZABOVE 

Ah 

Altitude  of  missile  above  sea  level  at 
t = 0 

RTD 

Conversion  factor  for  radians  to  degrees 

DTR 

Conversion  factor  for  degrees  to  radians 

RLBK 

L/B 

L/B,  actual  plume  ratio  of  length  to 
breadth 

NT 

Length  of  the  trajectory  time  and 
velocity  tables  (currently  NT  = 27) 

TMAS(K) 

Adjusted  target  time,  see  definition  of 

T 

* Data  input  from  card  reader 
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PROGRAM  SYMBOL 

VARIABLE 

DESCRIPTION 

199 

Flag,  if  the  current  launch  set  number 
is  a multiple  of  five  or  if  this  is 
the  first  launch  set  then  199  = 1. 
Otherwise  199  = 0 

188 

Flag,  if  this  is  last  launch  set  in 
automatic  selection  sequence  188  = 1, 
otherwise  188  = 0 

XE(I) 

*TE 

YE (I) 

V 

Target  acceleration  tables,  derived 

TE 

from  target  velocity  tables. 

ZE(I) 

^TE 

Z 

z^CO),  h 

Specific  launch  altitude  (from  launch 
set  table  ZTEO 

ZZZ,  ZALT 

^asL 

Specific  target  altitude  in  earth 
fixed  coordinates  + altitude  of 
missile  launch  site  above  sea  level 
(FT) 

AS 

zo 

Specific  launch  crossing  angle  (from 
launch  set  table  ASC(I) 

. 

R 

R. 

Specific  launch  initial  range  (from 

1 

launch  set  table  RR(I) 

XD 

Xte(°) 

YD 

Y^CO) 

Initial  target  velocities  in  earth 

fixed  coordinate  system 

ZD 

W0) 

ASCO 

*L 

Angle  from  positive  - axis  to  the 

vector  from  origin  to  target  (note 
this  is  just  -W) 

SCO 

*L 

in  DEG 

SO 

Initial  elevation  angle  of  target  (RAD) 

ASO 

0L 

0^  in  DEG 

2.16 
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PROGRAM  SYMBOL 

VARIABLE 

DESCRIPTION 

X 

X^g(O) 

X and  Y components  of  L in  earth 
^ os 

Y 

YteCO) 

fixed  coordinate  system 

SINSO 

sin(e^) 

COSSO 

cos(eL) 

SINSOO 

sin(ipL) 

COSSCO 

cos  0L) 

. . 

S2 

cos  (eL)cos  (i|/L) 

S3 

sin(eL)cos(4<L) 

S4 

cos(0L)sin(i);L) 

S5 

sin(0L)sin(^L) 

SPL 

sin(ipL) 

CTL 

cos(0L) 

CPL 

cos(l>L) 

STL 

sin(0L) 

KKK 

KKK  = 0,  if  Y^CO)  f 0 

1,  if  Y^CO)  = 0 

WW 

vT(0) 

Initial  velocity  of  target 

XDTGO 

x-tlCO) 

YDTGO 

^tlC0) 

Components  of  target  velocity 

in  launch  coordinates 

ZDTGO 

Ztl(0) 

SEO 

Dive  angle  (Note:  some  cases  are 
computed  in  SUBROUTINE  INIT) 
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III.  PROGRAM  MAIN  (Post-Real -Time) 


The  post -real -time  portion  of  the  STINGER  code  provides  an  estimate 
of  miss  distance,  lethality  and  simulation  quality.  In  addition,  the  post- 
real  -time  code  contains  portions  of  the  repeat  run  (rep-op)  logic. 

In  general,  the  analysis  in  post-real -time  data  will  contain  the 
following  events: 

• If  nonstationary  target,  calculate  Euler  angles  for 
transformation  from  launch-to-target  fixed  coordinates 

• Calculate  the  least  square  missile  approach  trajectory 
from  collected  A/D  data  and  project  miss  distance 
information 

• If  the  target  is  a high  speed  jet,  perform  a lethality 
test  with  the  probability-of-kill  model 

• Repeat  the  simulated  flight  20  times  and  obtain  average 
PK,  mean,  and  standard  deviation  of  miss  distance 
information 

• Give  control  to  pre-real  time  part  for  next  simulation 

More  information  concerning  post -real -time  data  processing  is 
contained  in  the  Section  titled  SUBROUTINE  MISCQM. 


TIME  CRITICAL  DATA  COLLECTION  TEST 


This  test  is  performed  to  verify  that  each  value  of  X^-X^collected 

in  real-time  corresponds  to  a value  in  the  range-to-go  tables.  The  test  is 
performed  by  testing  each  point  collected  as 

ML  < (Xj^-Xp^  < XXp  + DX,  i = 1,  NADJ 

where  )0C  is  the  i'th  entry  in  the  pre-camputed  range-to-go  table  and 
(XML_XIlPi  *s  actual  c^ata  P0*21*  collected. 


If  time  critical  data  collection  failure  is  detected  a diagnostic 
message  is  printed. 
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TABLE  3.1  DESCRIPTION  OF  MAIN  PROGRAM  SYMBOLS  (Post -Real -Time) 


i 

t 

I 


PROGRAM  SYMBOL 

VARIABLE 

DESCRIPTION 

NADJ 

The  actual  number  of  data  collection 
points  in  real-time 

XMISS(7) 

PPX(I) 

*ML_;XTL 

PPY(I) 

YML"YTL 

PPZ(I) 

ZML-ZTL 

TIME  (I) 

t 

Real  time  data  collection  points 

VMX(I) 

• • 

^ML~XTL 

VMY(I) 

YML'YTL 

VMZ(I) 

ZML‘ZTL 

XX(I) 

Range -data  collection  table 
generated  by  SUBROUTINE  INIT 

LEVEL 

Status  of  maneuver,  -7  implies  not 
in  real-time,  0 implies  in  real- 
time, +7  implies  target  trajectory 
table  exceeded 

ITT 

Target  trajectory  number,  note: 

ITT  = 1 is  a stationary  target 

ww 

VT(0) 

Initial  target  velocity 

IS1 

Print  flag,  IS1  f 1 indicates  print 
is  desired 

XMAN(1,I) 

t 

Target  maneuver  velocity  and  time 

• 

tables  at  range  table  data  collection 

XMAN(2,I) 

^TL 

point  (See  MAN(J,I)  in  SUBROUTINE 

• 

REALT) 

XMAN(3,I) 

ytl 

XMAN(4,I) 

^TL 

PROGRAM  SYMBOL 


VARIABLE 


DESCRIPTION 


XM(1)  ,YM(1),ZM(1) 
XM(2),YM(2),ZM(2) 
XM(3)  ,YM(3)  ,ZM(3) 
RM(1) 

RM(2) 

RM(3) 

AXM(ICN) 

AYM(ICN) 

AZM(ICN) 

ARM(ICN) 

AXB(ICN) 

AYB(ICN) 

AZB(ICN) 

ARB(ICN) 

AXC(ICN) 

AYC(ICN) 

AZC(ICN) 

ARC(ICN) 


Missile  trajectory  peircing  point 
in  Yp  - Zp  plane 

Missile  trajectory  piercing  point 
in  X£  - Y£  plane 

Missile  trajectory  piercing  point 
in  - Z£  plane 

Nearest  approach  of  missile 
trajectory  to  origin  in  the 
Yp  - Zp  plane 

Neaiest  approach  of  missile 
trajectory  to  origin  in  the 
X£  - Y ' plane 

Nearest  approach  of  missile 
trajectory  to  origin  in  the 
X£  - Z'  plane 


Storage  arrays  for  ICN'th  repeat 
run,  where  ICN  varys  from  1 to 
NRUNS.  The  computer  variables 
correspond  to  XM(1),  YM(1),  ZM(1) 
and  RM(1)  respectively 


Storage  arrays  for  ICN'th  repeat 
run,  where  ICN  varys  from  1 to 
NRLJNS.  The  computer  variables 
correspond  to  XM(2) , YM(2) , ZM(2) 
and  RM(2)  respectively 


Storage  arrays  for  ICN'th  repeat 
run,  where  ICN  varys  from  1 to 
NRUNS.  The  computer  variables 
correspond  to  XM(3)  , YM(3)  , ZM(3) 
and  RM(3)  respectively 


PROGRAM  SYMBOL 

VARIABLE 

DESCRIPTION 

YA 

See  equations  in  SUBROUTINE  MISCQM 

ZA 

a2 

See  equations  in  SUBROUTINE  MISCQM 

IV.  SUBROUTINE  FLIGHT  AND  REALT  (Real-Time) 


This  section  contains  a description  of  Subroutines  required  during 
real-time  operation  of  the  MICOM  Hybrid  STINGER  simulation.  The  documentation 
is  provided  in  three  parts: 


• Flow  charts 

• Equations 

• Description  of  Symbols 


The  flow  chart  given  in  Figure  4.1  only  to  the  digital  program. 
However,  the  equations  and  description  of  symbols  are  given  for  both  analog 
and  digital  implementations . 
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FLOW  CHART  FOR  SUBROUTINE  FLIGHT 
FIGURE  4.1 


M 


__  


sense  line  5 high 
indicates  static 
test  complete  on 
analog 


1 


A 


FLOW  CHART  FOR  SUBROUTINE  REALT 
FIGURE  4,2 
4.3 


t 


HAS 


YES 


f 

STORE  TIME  AND  TARGET 
VELOCITY  IN  A TABLE 


YES 
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PLUME  SHAPE 


7y  -Y  Tr 

K ML  TL; 


Tz  J7 

1 ML  TL; 


(XML-XTL^  + 


2 tan 


(DEG) 


(f)  00 

P_ INCH  + = _ (IN) 

2 | L | 

1 os 1 

P.RAD  + = tan'1  U/|L  j (RAD) 

4 Ub 

P2RAD+  = i/|Losl 


(RAD) 


PLUME  AZIMUTH  AND  ELEVATION 


*„  ■ * - P COS  (T)  (DEG) 

7 4 2 RP 

i>7  = (DEG/SEC) 

e7  = e4-P2  sin  (TRP)  (DEG) 


(DEG/SEC) 


PLUME  LENGTH  TO  BREADTH  RATIO  AND  ASPECT  ANGLE 


( ! 


( ) 


l ftj  L sin  (e)  if  B<L / B-3.0 


= L / 1-cos2 (e) 


l/B  = t / 1-cos2  (e) 
B ' 


cos  ( e ) = 


L • CT 
os  L 

L I | CT  | 
os  1 1 L 1 
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PLUME  IRIS  RATIO 


H 


c7  (7OOO  ) 


-2.341 


(H 


R=7000 

e 


(W/CM2 ) 


= P2INCH  SIN(Pl) 
2 COS^  (P./2) 


Jtu7=  K2At73,  K2  = .1,  K3  = 1.003258 


(IN2) 


(W/STER) 


I ) 


r7 


H _f2 
c7 

T 

f , plume 
Jtu 


From  IRSS  Calibration 
curve  fit 


17  17  (Ir7  } 


From  IRSS  Calibration  limits 


l 


r 


PLUME  ROTATION  ANGLE 


( ^ 


6 = cos 


-1 


/ L*  + Lz 
x y 


C = -cos  g 
x 


C = 

y 


+sin  g if  L > 0 
x 


-sin  g if  L <0 
x 


RP 


LZ  LZ 

CLY 

LY  "LX  * 

CLZ 

’LY  "LY  * 

CLX 

os  * V 

* C 

■ * cT  | 

os  L 1 

‘ 1' 

> 

! -e 

if  D 

: +e 

if  D 

• 3ir/2 

if  D 

-1 , 


C 


V Dz  + D2  + 

X Y Z 
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OS  x CL 


FIGURE  4.6-i- 


APPARENT  TARGET  PLUME  AS  SEEN  FROM  THE  MISSILE  SEEKER 
(THE  PLANE  OF  THE  PAPER  IS  PERPENDICULAR  TO  THE  L VECTOR) 


TARGET  ANGLE  OF  ATTACK  AND  MACH  NUMBER 


- . 00003h 


p = .00237692  e 


(SLUGS/FT3) 


VT  ^ *TE  + *TE  + ZTE 


(FT/SEC) 


C = CLAA 

La  IQm 


a = 1116.89  - .003894  h _T  (FT/SEC) 

S clSL 


Q = V /a 
M T S 


/ X2  + Z2  sin  ("  tan-1 
ML  ML  1_ 


~ZML 

) 

ML 


+ \] 


+ Ah*  (FT) 


— Altitude  of  launch  site  above  sea  level 


pV|  CLSTg 
a 

T 


/ X2 
TE 


+ Y2  + Z2  ) 
TE  TE  / 


(RAD) 


*THIS  MODEL  IS  NOT  CURRENTLY  INCORPORATED  IN  THE  MICOM  REAL 
TIME  SIMULATIONS. 


FLARE  TRANSPARENCY  RADIUS  AND  IRIS  RATIO 


Jtu2  = Jtu2(tFD>  ' (W/STER),  DFG 


I ) 


RF  " / (XML-X2L)2  + (V-X2L)2  + <ZML-Z2L)2 


Hc2  _ Jtu2  ' RF  (W/CM2 ) , DFG 


fc2  “ *2  (HC2» 


*2  " *2  «c2> 


(RAD) , DFG 


(RAD)  , DFG  based  on  IRSS  calibration  limits 


J 


FLARE  AZIMUTH  AND  ELEVATION 


°Y2  " YML~Y2L 

X -X 
ML  2L 


(RAD) 


_ I 57. 3 (o  -o  ) , t>t 
J24REL  " Y4  Y2  “ EJECT 

°'  t<fcEJECT 


*^24REL  57*3|pYML  Y4LJ  °Y4  (XML  X4L)- 


X -X 
ML  4L 


r(YML~Y2L)  °Y2(XML~X2L) 
L XML-X2L 


*2  = V *24REL  (DEG) 


*2  = *4  -'"24REL  ^^/SEC) 


(DEG.) 


(DEG/SEC) 


-57 • 3 (° z4-° Z2 J ' t>tEJECT  (DEG) 


24REL 


0,  t<t 


EJECT 


624REL  "57*3 


(ZML  Z4L}  ~ °Z4(XML  X4L} 
• f\L-X4L) 


r(ZML  Z2LJ  ” °Z2(XML~X2L} 

L (X  -X  ) 

ML  2L 


(DEG/SEC) 


4\22 


FLARE  DYNAMICS 


( ) 


X2L  " 

-P 

CD2  S2V2 

i 

(N 

X 

g sin  eL  , t>tEJECT 

2M2 

Y2L  " 

“P 

CD2S2  V2 

X2L  ' 

t>t 

EJECT 

2M 

2 

Z2L  = 

-P 

CD2S2  V2 

Z2L  + 

g cos  e , t>t 

Li 

(FT/SEC2) 


2M, 


X2L=  X4L 


Y = Y 
2L  4L 


Z2L  Z4L 


X2L=  X4L 


Y = Y 
2L  *4L 


Z2L  Z4L 


t<t 


EJECT 


(FT) 


TAILPIPE  AZIMUTH  AND  ELEVATION 


°Y4  YML~Y4L 
XML_X4L 


(RAD) 


U = 57,3  °Y4 


(DEG) 


57.3  f (Y  -Y  ) - a (X  -X  )- 
| ML  4L  Y4  ML  4L 

(X  -X  ) 

ML  4L 


(DEG/SEC) 


ZML~Z4L 

XML-X4L 


(RAD) 


0 . = -57.3  q . 

4 Z4 


(DEG) 


6 = -57.3  r (Z  -Z  ) - o (X  -X  ) 

4 ML  4L Z4  ML  4L 

(XML-X4LJ 


(DEG/SEC) 


PITCH,  YAW  AND  ROLL  POSITIONS 


P = 57.29  / q'dt 


Y = 57.29  fx'dt 


(DEG) 


R = R(p') 


(DEG) 


4,24 


TABLE  4.1 


DESCRIPTION  OF  PROGRAM  SYMBOLS 
FOR  SUBROUTINE  FLIGIT  AND  REALT 


VARIABLE 

PROGRAM 

SYMBOL 

DESCRIPTION 

SPO 

Scaled  plume  rotation  angle 

trp 

TRP 

Plume  rotation  angle  (RAD.) 

6L 

THETAL 

Elevation  angle  of  target  (DEG.) 

t 

DT 

Real  time  (SEC.) 

g 

Gravitational  constant 

*E 

XE 

( Target  aerodynamic  accelera- 
1 tion  table  in  earth  fixed 

YE 

/ coordinates  (Note  ZE  includes 

£ 

] an  added  factor  of  32.174 

ZE 

ZE 

| FT/SEC/SEC 

a 

A1 

Target  angle  of  attack  (RAD.) 

f 

G , GAM 

Time  at  approximately  1000 
feet  to  go 

CL 

“t 

CLA 

Aerodynamic  lift  coefficient 
of  target  due  to  angle  of  attack 

V 

TI 

VTI 

Target  inertial  velocity 

l/b 

RLB 

Apparent  plume  length  to 
breadth  ratio 

E 

Subscript  which  denotes 

earth  fixed  coordinates 

*L 

Azimuth  angle  of  target  (DEG.) 

VARIABLE 


PROGRAM 

SYMBOL 


DESCRIPTION 


Subscript  which  denotes 
target  fixed  coordinates 


Subscript  which  denotes 
launch  coordinates 


Initial  range  (FT.) 


Range  (FT) 


Scaled  y/t'. 


Air  density  (SLUGS/FT.3) 

Minima  of  the  overload 
function 


Variable  scale  factor. 


K=  ( 1 + t 1 , t<  t' 


1 + Y , t>t’ 


Range  table  for  data  collection 


Scaled  range  table  (XX) 


28 


VARIABLE 

PROGRAM 

SYMBOL 

DESCRIPTION 

J 

. I 

r7 

RN 

For  MICOM  Hybrid,  a uniform 
random  number,  RN  e (-1,1) 

For  IRSS , iris  ratio  number 

7 

Cos  (e) 

COSE 

Cosine  of  angle  between  LOS 
and  center  line  of  target 

x -x 

ML  TL 

y -y 

ML  TL 

DX 

DY 

\ Scaled  XXX,  YYY,  ZZZ  (NOTE:  these 
) variables  may  overload  on  the 
( AD/4  until  time  = t'); 

z -z 

ML  TL 

DZ 

( 

t 

DT 

Real  time 

X -x 

ML  TL 

XXX 

x-Missile  minus  x-target 
position  (FT.) 

y -y 

ML  TL 

YYY 

y-Missile  minus  y-target 
position  (FT.) 

z -z 

ML  TL 

ZZZ 

z-Missile  minus  z-target 
position  (FT.) 

• • 

x -x 

ML  TL 

XDOT 

x-Missile  minus  x-target 
velocity  (FT. /SEC.) 

• • 

y -y 

ML  TL 

YDOT 

y-Missile  minus  y-target 
velocity  (FT. /SEC.) 

• • 

z -z 

ML  TL 

ZDOT 

z-Missile  minus  z-target 
velocity  (FT. /SEC.) 

4.29 


h ZALT  Altitude  of  missile 

above  sea  level  (FT.) 


XCOMP  Interpolated  value  of 

XDTGMS 


YCOMP  Interpolated  value  of 

| YDTGMS 


ZCOMP  I Interpolated  value  of 

I ZDTGMS 


XC  I Interpolated  value  of 

I XDM 


YC  I Interpolated  value  of 

| YDM 


ZC 


Interpolated  value  of 
ZDM 


PROGRAM 

SYMBOL 


DESCRIPTION 


XDM 

YDM 

ZDM 

TIME 

PPX 

PPY 

PPZ 

VMX 

VMY 

VMZ 


Tables, of  target  velocity 
components  in  earth  system 


Storage  for  missile 
position,  velocity  and 
time  in  the  region  of 
pre-specif ied  range 
table  entries 


IPTS 

MISSED 


Number  of  points  in 
range  table 

An  array  of  values 
representing  miss 
conditions 


Till 


Angle  between  (LOsxCL)and 
horizon  reference  vector 


VARIABLE 


PROGRAM 

SYMBOL 


x-Missile  minus  x-target 
(in  launch  coordinates) 


y-Missile  minus  y-target 
(in  launch  coordinates) 


z-Missile  minus  z-target 
(in  launch  coordinates) 


Target  inertial  velocity 


Interpolated  target 
accelerations  in  earth 
fixed  coordinate  system 


Longitudinal  center  line 
of  target 


Components  of  C in  earth 
L 

fixed  coordinates 


Line  of  sight  vector 


4.32 


PROGRAM 

SYMBOL 


DESCRIPTION 


FI 

F2 

F3 


Components  of  Lqs  in 
earth  fixed  coordinates 


x-Missile  position  in 
launch  coordinates 


y-Missile  position  in 
launch  coordinates 


z-Missile  position  in 
launch  coordinates 


x-target  position  in 
launch  coordinates 


y-target  position  in 
launch  coordinates 


z-target  position  in 
launch  coordinates 


x-position  of  target  in 
inertial  coordinates 


y-position  of  target  in 
inertial  coordinates 


z-position  of  target  in 
inertial  coordinates 


PROGRAM 

SYMBOL 

DESCRIPTION 

XDTGO 

i 

/ Initial  values  of  XDTGMS, 

I YDTGMS , ZDTGMS 

YDTGO 

J 

ZDTGO 

Commanded  irradiance  for 
plume 

19.5 

Focal  length  of  projection 
lens  (CM) 

HR7E 

Irradiance  at  7000  meters 
(W/CM2) 

PlIRSS 

Vertex  angle  of  plume 
transparency  (DEG) 

P2IRSS 

Length  of  plume  trans- 
parency (IN.) 

AT  7 

Area  of  plume  trans- 
parency (IN?) 

PJTU1 

Available  radiant 
intensity  of  target  as 
function  of  plume  trans- 
parency area  taken  from 
calibration  in  IRSS  (W/ 
STER) 

0.3 


Neutral  density  trans- 
mission factor  for  plume 


VARIABLE 

PROGRAM 

SYMBOL 

DESCRIPTION 

Wm 

T 

WT 

Weight  of  target 

Sm 

ST 

i 

Reference  surface  area  of 

T 

target 

c 

Horizon  reference  vector 

°M 

QM 

Mach  number  of  target 

a 

s 

AAA 

Speed  of  sound 

V2 

Speed  of  flare 

R 

F.ange  to  go  between  missile 

F 

and  flare 

H.  ,H 

Irradiance  of  flare  W/CM2 

i2  c2 

(Note,  letter  c denotes 

calibration) 

• 

X 

/ Velocity  of  flare  in 

FL 

I launch  coordinates  (FT/SEC) 

• 

y 

) 

FL 

) 

• 

Z 

( 

FL 

' 

°Y4 

Approximate  ^ (RAD) 

°Z4 

Approximate  (RAD) 

°Y2 

o 

YT 

Approximate  (RAD) 

4-35 


PROGRAM 

SYMBOL 


DESCRIPTION 


°ZT 


J 


tu2 


fcFD ' fcf d 


*^F2 

*T4 

*T7 

it 

T2REL 


Approximate  0 (RAD) 

i 

Available  radiant  intensity 
of  flare 

Time  of  flare  drop  (SEG) 

Norm  of  line  of  sight  (FT) 

Time  associated  with  flare, 
drop 

Missile  wing  deflection 
Flare  azimuth  (DEG) 

Tail  pipe  azimuth  (DEG) 
Plume  azimuth  (DEG) 

Flare  elevation  (DEG) 

Tail  pipe  elevation  (DEG) 
Plume  elevation  (DEG) 


Flare  azimuth  rate  (DEG/SEC) 

Tail  pipe  azimuth  rate  (DEG/ 

SEC) 


VARIABLE 


PROGRAM 

SYMBOL 


DESCRIPTION 


*7 

*T7 

Plume  azimuth  rate  (DEG/SEC) 

*24 

• 

*T2REL 

l 

' 

• 

62 

®T2 

Flare  elevation  rate  (DEG/SEC) 

e. 

®T4 

Tail  pipe  elevation  rate 

4 

(DEG/SEC) 

°7 

°T7 

Plume  elevation  rate  (DEG/SEC) 

*24 

• 

6T2REL 

fc2 

Transparency  radius 
for  flare 

used  to 
command 

fc4 

Transparency  radius 
for  tailpipe 

true 

size 

P 

6 

GUM  pitch  angle  (DEG) 

Y 

. 

* 

GUM  yaw  angle  (DEG) 

R 

* 

GUM  roll  angle  (DEG) 

r’ 

r 

GUM  Yaw  rate  (DEG/SEC) 

q' 

GUM  pitch  rate  (DEG/SEC) 

p' 

GUM  roll  rate  (DEG/SEC) 

4.37 


PROGRAM 

SYMBOL 

NPATH 

LOOP1 

LOOP2 

INFLRT 

LTF 

NSEG 

IT1 

INDEX 

LEVEL 

IT1,  IT2 


Flag  to  control  entrance  into 
computation  loop  of  real-time 
code. 


Loop  counter  for  real-time 
interrupts  prior  to  sense  line 
6 coming  high 

Loop  counter  for  real-time 
interrupts  after  sense  line 
6 comes  high 

Flag  to  signal  entrance  into  real 
time,  INFLRT  = -1  indicates  batch 
and  INFLRT  = +1  indicates  entrance 
into  real-time 

Clock  time  at  beginning  of  real- 
time loop.  This  value  is  used 
, later  to  compute  specific  loop 
time  for  each  loop 

Segment  number  of  real-time  code 
which  is  to  be  entered. 


Index  for  storing  real-time  data 
for  post  real-time  analysis 

Status  of  maneuver 
-7  -*•  not  in  real-time 
0 in  real-time 
+7  -»■  target  trajectory  table 
exceeded 

Index  used  in  interpolation  of 
target  trajectory  time 


NT 


Length  for  target  trajectory  tables 


PROGRAM 

SYMBOL 

DESCRIPTION 

DXG 

■ 

DYG 

DZG 

Dummy  storage  location  for  current 
value  of  DX,  DY,  DZ,  XDOT,  YDOT 

XDO 

and  ZDOT  respectively.  These 
variables  are  returned  to  the  main 

YDO 

program  for  diagnostic  purposes 

ZDO 

PL 

Actual  plume  length  (FT) 

PB 

Actual  plume  breadth  at  base  (FT) 

PLA 

Apparent  plume  length 

t 

MAN (1,1) 

^TL 

MAN(2,I) 

Target  time  and 

ytl 

MAN (3, I)  velocity  at  each 

range  table  data 

ZTL 

MAN(4,I)  collection  point 

1 

V.  REAL-TIME  I/O 


This  Section  contains  a description  of  real-time  input  output 
requirements  of  the  MI  COM  Hybrid  STINGER  Simulation.  Hie  tables  and  ; 

figures  included  in  this  section  describe  the  flow  of  analog  and  discrete 
information  between  an  AD/4  analog  computer,  a CDC/6600  digital  computer 
and  the  IRSS. 


I - 

NOTE:  The  CTRL  and  SL  word  are  also  used  to  transmit  SHDCS 
and  DIs,  the  specific  usage  is: 


o 

SHDC 

8 = 

CTRL  8 

o 

DI 

3 = SL  3 

o 

SHDC 

9 = 

CTRL  9 

/ o 

DI 

7 = SL  7 

o 

SHDC 

10 

= SL  11 

o 

DI 

8 = SL  8 

o 

SHDC 

11 

= CTRL  11 

o 

DI 

10  = SL  10 

FIGURE  5.1 


DISCRETE  CIRCULATION 


TABLE  5 . 1 

DISCRETE  SIGNALS  FROM  AD/4  TO  CDC/6600 


AD/4 

PATCH 


CDC/6600 

BIT 


ORIGIN  ULTIMATE 
DESTINA- 
TION 


SL11  NORMAL  TERMINATION  (SHDCIO) , SIGNALS  RAMP  AD/4  IRSS/19 

DOWN  COMPLETE 


TR27  SLO  LOOP  TO  BATCH,  RETRANSMITTS  ICs  WHEN  SET  AD/4  6600 

LOW 

TR22  SL5  STATIC  TEST  OK,  SIGNALS  AD/4  READY  TO  AD/4  6600 

RECEIVE  IC 


RAMP  UP  READY,  SET  HIGH  BY  OPERATE  SWITCH  AD/4 
IN  MICOM  HYBRID,  


THIS  STARTS  REAL  TIME  LOOPING 

6600  ABNORMAL  ABORT,  USED  TO  EXIT  REAL  AD/4 
TIME  LOOPING 


IRSS  READY  (DI3) 

ND  OF  PROBLEM  (DI7) 


RGENCY  SHUTDOWN  (DI8) 


INITIAL  FIELD  REQUEST  (DUO) 


TABLE  5.2  DISCRETE  SIGNALS  FROM  AD/4  TO  IRSS 


AD/4 

PATCH 

IRSS 

BIT 

DESCRIPTION 

ORIGIN 

ULTIMATE 

DESTINA- 

TION 

T46 

SHDCO 

SHUTTER  #2  COMMAND 

AD/4 

IRSS/3 

T61 

SHDC1 

SHUTTER  #4  COMMAND 

AD/4 

IRSS/3 

T63 

SHDC2 

SHUTTER  #7  COMMAND 

AD/4 

IRSS/3 

SHDC8 

SIMULATION  RUNNING,  (THIS  PATH  NOT  USED 

SEE  CTRL 8) 

SHDC9 

FIELD  TRANSMITTED,  (THIS  PATH  NOT  USED 

SEE  CTRL9) 

SHDC10 

NORMAL  TERMINATION,  (THIS  PATH  NOT  USED' 

SEE  SL11) 

T43 

SHDC12 

TARGET  AQUIRED 

AD/4 

IRSS/3 

T40 

SHDC13 

DISCRETE  #5,  SWITCH  SEEKER  FROM  EXTERNAL 

AD/4 

IRSS/3 

> 

TO  INTERNAL  POWER 

NOTE:  . IRSS/3  IS  IRSS  CONSOLE  3 AND  IRSS/19  IS  IRSS  CONSOLE  19 f THE  DIRECT 
CELL  DATA  LINK 


TABLE  5.3  DISCRETE  SIGNALS  FROM  IRSS  TO  AD/4 


AD/4 

PATCH 

IRSS  BIT 

DESCRIPTION 

ORIGIN 

ULTIMATE 

DESTINA- 

TION 

T72 

DI3 

IRSS  READY,  NOT  SENT  UNTIL  TARGET  ACQUIRED 
SIGNAL  IS  RECEIVED 

IRSS 

6600 

T56 

DI7 

END  OF  PROBLEM,  SIGNALS  EQUIPMENT  SHUT 

DOWN  IN  IRSS 

IRSS 

AD/4  660( 

T55 

T53 

DI8 

DI10 

EMERGENCY  SHUT  DOWN,  SIGNALS  EXCESSIVE 
ACCELERATION  IN  GUM 

INITIAL  FIELD  REQUEST 

IRSS 

IRSS 

AD/4  660( 

6600 

T53 


110 


[INITIAL  FIELD  REQUEST 


IRSS 


6600 


DISCRETE  SIGNALS  FROM  CDC/6600  TO  AD/4  OR  IRSS 


AD/4  I CDC/660 

PATCH  I B 


DESCRIPTION 


ORIGIN 


TR33  I CTRL4  | ICs  SENT 


CTRL7 


CTRL15 


• 

1 

INTECEPT  DETECTED,  SET  HIGH  BY  INTERCEPT 

OR  MISS  CONDITION 

rz 

1 

• 

REAL  TIME  LOOPING,  SIGNAL! 

3 CDC/6600  WAIT- 

ING  FOR  SL6  HIGH 


CTRL9  FIELD  TRANSMITTED  (SHDC9)  SIGNAL  TO  GEPAC  6600 
THAT  FIRST  FRAME  OF  DATA  IS  AVAILABLE. 

THIS  OCCURS  A* TER  THE  6600  HAS  RECEIVED 
INITIAL  FIELD  REQUEST,  DUO 

CTRL 8 SIMULATION  RUNNING  (SHDC8)  TURNS  ON  INDI-  6600 

CATOR  LIGHTS  ONLY.  SENT  WITH  SHDC11 


CTRL11  START  SIMULATION  (SHDC11)  IN  CLOSED-LOOP  66C0 
GEPAC  RESPONDS  BY  SENDING  FIRST  CLOCK 
PULSE.  SECT  WITH  SHDC8 


5 


TARGET 

GENERATION 


FIGURE  5.3  ANALOG  SIGNAL  CIRCULATION 


tf'DLuiuJ.  t^/u.v  ^ - acifiDla 


- ADCS-  to  QEPi 


i 


(ADC  ASSIGNMENT  CONTINUED) 


5.10 


zv  s 


TABLE  5.8  COMMANDS  TO  IRSS/9 
V17  - U02+-W9X28A 


ANALOG  VARIABLE 

NAME 

AD/4  ANALOG 
VOLTAGE 

RELATIVE  TRUNK 

ASSIGNMENTS 

AD/4  TRUNK  LINE 
ASSIGNMENTS 

r' (YAW  RATE) 

1.0  r’ 

5 

354 

q’  (PITCH  RATE) 

1.0  q' 

6 

355 

p’  (ROLL  RATE) 

-3 

13.8x10  p' 

7 

356 

1 

trp 

.5555  T^p 

s 

2 

351 

k 


AD/4  TO  C DC/6600 

(sent  to  IRSS  via  direct  cell) 

V12-*- 


AD/4  ANALOG 
VOLTAGE 

RELATIVE  ADC 

ASSIGNMENT 

AD/4  TRUNK  LINE 
ASSIGNMENT 

S.53X103t2 

1 

110 

9.534X103t4 

2 

111 

9.896X102P2 

3 

112 

.5555  R 

'4' 

113 

1.7255 

5 

114 

10.742  i2 

6 

115 

10.742  i4 

7 

116 

10.741  i? 

8 

117 

1.1111  ii2 

9 

130 

3.333  02 

10 

131 

1.111 

2 

11 

y 

132 

3.333  04 

12 

133 

1.111  \l> 

7 

13 

134 

3.333  6 

7 

14 

135 

1.25  P 

15 

136 

1.111  Y 


16 


137 


i 


DYNAMICS.  INC. 


VI.  SUBROUTINE  KSCALE 


This  subroutine  estimates  t',  the  time  at  1000  feet-to-go  for  target 
missile  intercept  and  computes  the  overload  function  for  real-time  variable 
scale  factors.  The  estimates  are  based  on  a simplified  model  involving  a 
constant  bearing  course  and  a nomial  missile/ target  trajectory. 

In  addition,  the  program  computes  y min  - the  minimal  of  the 
overload  function,  t y min  - the  time  at  which  y min  occurs  and  the  smallest 
bit  ADC  resolution  for  relative  missile  positions. 


EQUATIONS: 


• Qg  = tan 


•1  / »'  rTL(0)2+ZTL(0>' 

R.-1000+t'x  (0) 

1 1 L 


• R.-A  R cos  ( eo)  - 1000  + t'  x T (0)  = 0 

1 Z 1 .Lf 


(Note,  the  interative  solution  of  this  equation  yields  t') 


or 


• K = -xT  , /R.  = +20,  000/R.  -*•  K R.  = 20,  000 

p Llm/  i i p i 

fmin  ty  min 


• K R.  =(  1+ 
P 


■)( 


,+x  (0)  t - AR  . cos  6, 


TL  ymm 


min 


y min 


20,  000  t' 


t y min 


(R +x_T  (0)  t . - AR  . cos  e ) 

i TL  ymin  mm  2 


- t' 


dy 


min 


= At  . + Bt  . + C = 0 
y min  ymm 

y min 

(Note,  the  interative  solution  of  this  equation  yields  t 


y min^ 


• a.*tl<0) 


• B = 2 x_T  (0)  ( 20,  000  - R + AR  . cos  0O)  - 20,  000  V cos  ( 
1 J_j  min  z m 


C = (R  - Ar  . cos  6 _)  (20,000  - R + Ar  . Cos  e_) 
min  2 min  2 


'min 


- (v* 


20, 000 


X (0)  t y . - A R cos 
TL  ' min 


- \( 

e2  +1  M t Ymin  ) 


A2  _ (A Voltage)  (Analog  Scale  Factor  for  DX) 

1 + Y. 


min 


Here  A Voltage  is  the  analog  voltage  change  resulting  from  a smallest 
bit  change  in  the  associated  DAC. 


A V = .024 

Analog  Scale  Factor  = 20 


Typical 

Values 


■ 


Subroutine 

KSCALE 


FIGURE  6.3 

OVERVIEW  OF  SUBROUTINE  KSCALE 


6.4 


RANGE  VOLTAGE  7 - SELECTION 


SELECTION  CRITERIA 


i 


FIGURE  6.6 

CONVERGING  SEARCH  FOR  y . 


6,7 


500 


20  TIME  (SEC) 


VELOCITY  VS  TIME  FOR 
NOMINAL  STINGER  TRAJECTORY 

FIGURE  6.7 


6.8 


r 


PROGRAM  VARIABLE 


SYMBOL 


DESCRIPTION 


EDOT 

Y/t' 

Y/t' 

RI 

Ri 

Initial  range  (FT) 

G 

t' 

Time  at  approximately  1000  feet  to 
go 

GG 

^Yinin 

GGG 

Ymin 

Minima  of  the  overload  function. 

obtained  by  iterative  solution 

IS1 

Print  flag,  IS1  j1  1 indicates  print 
is  desired 

THETA1 

el 

See  Figure  titled  Simplified 

Geometry  for  Estimating  Target- 

THETA2 

02 

Missile  intercept. 

XDTGO 

XjL (0) 

Initial  components  of  target 

• 

YDTGO 

Ytl(O) 

velocity  in  launch  coordinates 

ZDTGO 

z^CO) 

DG 

At’>AtYmin 

Time  increment  in  t*  search 

AN 

Scale  factor  for  decrementing  At 
as  the  search  is  refined 

DELTAR(I) 

“(V 

Distance  missile  has  flown  along 
nominal  trajectory 

VM(I) 

VV 

Velocity  of  missile  for  a nominal 
trajectory 

TAMA(I) 

Sn 

Time  table  for  the  above  nominal 
trajectories 

“p 

Maximum  analog  voltage  = lOOv 

NB 

Length  of  nominal  trajectory  tables 
TAMA,  VM  and  DELTAR 

XLlm 

Maximun  value  expected  for  (used 

to  estimate  scale  factors,  XVTw  = 
20,000  FT) 

Maximum  analog  voltage  (k^  = lOOv) 

Smallest  bit  ADC  resolution  in 
feet 


VII.  SUBROUTINE  MISCOM 


DYNAMICS.  INC 


VII.  SUBROUTINE  MISCOM 


Subroutine  MISCOM  transforms  the  relative  missile  minus  target 
location  in  launch  coordinates  to  target  fixed  coordinates  referenced  to 
the  target  tailpipe.  The  missile  path  projection  is  then  computed  by  a 
least  square  straight  line  approximation  of  real  time  data.  Nearest 
approach  miss  distance  is  computed  in  each  plane  of  the  target  fixed 
coordinate  system  (referenced  to  tailpipe  origin). 


H 


I 

■ 


MISSILE  PATH  PROJECTION 


The  missile  path  projection  is  computed  in  target  fixed  coordinates 
referenced  to  the  target  tailpipe.  This  projection  is  based  upon  a linear 
curve  fit  of  real  time  relative-position -data. 

For  a least  squares  approximate  linear  curve  fit,  the  projection 
equations  are 


Yp  = a2  Xp  +b1 


Z”  = a2  +b2 


The  least  square  approximate  fit  for  a^  and  b^  is  obtained  from  the  simultaneous 
solution  of  the  following  equations. 

n n n 

+ bi.£,V  -,£  V V 

1=1  1 1=1  1 1=1  1 1 


n 


n 


a E X + b n = Z Y_ 
Xi=l  Fi  1 i=l  Fi 

Then,  and  f>2  can  be  solved  similarly 


n 2 

n 

+ 

b2.ZXF 

1=1  l 

i=l  i 

n 

n 

+ 

bn  = E z 

1=1  1 

i=l 

where  n is  the  number  of  real  time  data  collection  points. 


Then,  the  piercing  points  of  the  missile  trajectory  into  the  - Y ' plane,  X^  - 
ZF  plane  and  the  YF  - ZF  plane  is  computed.  These  values  are  printed  as: 


XM  (1) , YM  (1) , ZM  (1)  - the  piercing  point  in  the  Y^,  - Z^,  plane 


Then,  the  nearest  approach  of  the  missile  trajectory  to  origin  is  computed  in 
th®  Yp  - Zp  plane,  Xp-  Yp  plane  and  the  Xp  - Zf  plane  respectively  as: 


RM(1)  ~ V XI2  + Y'2  + Z'2 
yi  *fi  + zfi 


FM  (2)  — J y 1 2 x y 1 2 7 » 1 2 

XF2  + YF2  + v2 


M(3>  ’ ''K’? 

F3  F3  F3 


Here  the  origin  of  the  X^,  - Y^,  - Z^,  system  is 


“ - 

— « 

X' 

0 

FI 

Y* 

0 

F3 

Z’ 

0 

F2 

m m 

■ — 

TABLE  7.1  DESCRIPTION  OF  VARIABLES  IN  SUBROUTINE  MISCQM 


S 


COMPUTER 

VARIABLE 


PHYSICAL  VARIABLE 


DESCRIPTION 


S2 

cos  (0L)  cos 

Elements  of  coordinate  transformation 
launch  coordinates  ±o  earth  fixed 

SPL 

coordinates 

S3 

sin  (0L)  cos  (ipL) 

S4 

cos  (0L)  sin  (i^L) 

CPL 

cos  (i|jl) 

• 

S5 

sin(0L)sin(i()L) 

STL 

sin(0L) 

CTL 

cos (0L) 

N 

The  number  of  data  points  collected 
during  real  time  data  collection 

WW 

vT(0) 

Speed  of  target  at  time  = 0 

XM(1)  ,YM(1)  , 

Y • V • *7  • 

aF1 ' *F1 ' “FI 

Missile  trajectory  piercing  point  in 

ZM  (1) 

YF  " zf  plane 

XM  (2)  , YM  (2)  , 

XF2'YF2'ZF2 

Missile  trajectory  piercing  point  in 

ZM  (2) 

1 

Xp  ~ Yp  plane 

XM  (3)  ,YM  (3) 

XF3'YF3,ZF3 

Missile  trajectory  piercing  point  in 

ZM  (3) 

RM(1)  , RM ( 2 ) , 

Xp  - Zp  plane 

Distance  of  each  piercing  point  from 

RM  (3) 

origin  of  X^  - Y^  - Z^  system 

YA 

al 

See  equations 

ZA 

a2 

See  equations 

VIII.  SUBROUTINE  GUID 


(Coordinate  transformation  from  earth  fixed  to  launch  coordinates) 


X “ 

coseL  cos^ 

coseL  sin^ 

-sine^  ' 

X' 

yl 

= 

-sijnfiL 

cos^L 

0 

ye 

ZL 

sineL  cos^ 

sine^  sim|^ 

COS0L 

ZE 

■ 

8.0 


TABLE 8.1  DESCRIPTION  OF  SYMBOLS  FOR  SUBROUTINE  GUID 


A 


P 


PROGRAM  VARIABLE 

SYMBOL 

DESCRIPTION 

i 

XD 

Components  in  earth  fixed  coordinate 

YD 

system 

ZD 

XDTG 

Components  in  launch  coordinate 

YDTG 

system 

ZDTG 

■ 

OOSSO 

+ coseL 

cossco 

+ cosily 

SINSCO 

+ sim|>L 

SINSO 

+ sine^ 

r 


— 


1 


IX.  SUBROUTINE  INIT 


This  subroutine  computes  the  range  table  for  real-time  data 
collection  and  the  initial  conditions  for  the  analog  computer.  More 
* detail  concerning  the  equations,  given  here,  can  be  obtained  in  the 

section  dealing  with  real-time  digital  computation. 


Initial  Velocity  of  Target  in  Launch  Coordinate  System 


m m 

coseL  cosfL  COS0L  "SineL 

Xte(0) 

Yu(0) 

-sin^  cos^  0 

>(0>  . 

sine^  cos^L  sineL  coseL 

m m 

a ■ 

/ — \ 

O 
' — / 

where  e = eT (0) 

L Li 

*L  = *L(0) 

This  transformation  is  performed  in  SUBROUTINE  GUID 

NOTE:  The  transformation  is  orthonormal  therefore  A'^  = AT 

Dive  Angle 

Zrrrc(O) 

SEO  = tan'1  - . , if  Y^CO)  t 0 

^(0)  + Y^(0) 

NOTE:  For  all  other  cases  see  computations  performed  in  program  MAIN 

Initial  Angle  of  Attack 
a = 1116.89  + .003894  h (SLUGS/FT3) 

oSL 

p = .00237692  e-00003hasL  (FT/SEC) 

Qm  " Vas 


- CLM 


a = 


2Wt 


pvf  cL 

°T 


g 


x|e(0)  + %(0)  + %(0) 


Initial  cosfe) 

cos(«)  = t(sineL)  (injCO)  cos«  - ^(0)  * * (0)  s“a) 

- (cosa  . ZteCOJ  smc/  ^(0)  ♦ y^toj)  • 

• • 

(XpgCO)  coseL  cosily  + Y^fO)  coseL  sin<(>L)}/VT 

Initial  Apparent  Plume  Length  to  Breadth  Ratio 
&/B  = (L/B)  ^l  - cos2(e) 


Initial  L 

— ~ os 

1^  = -R^  cos6l  cosii>L 
Ly  = -R^  cos0L  sim|>L 
Lz  = R sin0L 


9.1 


Initial  Horizon  Reference  Vector 


L 

e = cos'-*-  ■ 

*^2  + L2 
x y 


Cx  = -cos  (3) 

Cy  = -sin  (3),  if  XLqs  < 0 
+sin  (3) , if  XLqs  > 0 


Initial  Target  Centerline  Vector 


= Xj^CO)  fcosa  + Z^fO)  sino) 


/J2 


c,  - Y^fO)  (coso  * ZT£(0)  sine) 
/•Z 


Xfe(O)  + Yfe(O) 


CL  = Z^CO)  COSO  - / 


X^fO)  * Y^CO) 


sma 


Initial  Plume  Rotation  Angle 
6 -cos-1  Cx  (LyClz  • LzCy)  + Cy(LA  - LA) 

Vlz  ’ LA/*  ^ - W2*  - W2 

Dz  = LxCLy'LyCLx 


9,2 


* 


r = ir/2  - e , if  D < 0 
rp  z 


tt / 2 + 0 , if  D >0  and  e < w/2 

Z 


e-3ir/2,  if  D > 0 and  0 > tt/2 
z 


Initial  Plume  Iris  Ratio 


H = 1209.675  H 
c/ 


R=7000 


- 2 . 341 (ln(R./22965 . 831) ) 
“ 1 


P,  - 2tan'W) 


P - (f)  CO 


m»i 


At7  " p2  sln  ffp 
2cos2(P1/2) 


T _ 1. 003258 (lnA.,/10) 
Jtu7  " e 


H 


c7 


r7 


lfplume 

Jtu 


X1  W 


(w/cm2) 

(DEG) 

(IN) 


(IN)' 


(W/STER) 


(table  look  up) 


9.3 


TABLE  9.1  DESCRIPTION  OF  SYMBOLS  FOR  SUBROUTINE  INIT 


PROGRAM  VARIABLE 


SYMBOL 


DESCRIPTION 


Number  of  real-time  data  collection 
points 

Distance  between  target  and  missile 
at  last  data  collection  point 

Distance  between  real-time  data 
collection  points 

Minima  of  the  overload  function 
(see  SUBROUTINE  KSCALE) 

The  difference  between  the  data 
collection  distance  commanded  and 
the  largest  data  collection  distance 
possible  with  the  analog  scaling 
given. 

New  value  for  NPTS  if  NPTS  is  too 
large 


Target  initial  velocity  in  earth 
coordinates 


Target  initial  velocity  in  launch 
coordinates 


KKK  = 0 if  Y^CO)  t 0 
1 if  Y^CO)  = 0 

Dive  angle  (Note,  some  cases  are 
computed  in  MAIN) 


PROGRAM  VARIABLE 

SYMBOL 

DESCRIPTION 

XDTGM(I) 

YDTGM(I) 

XTL 

^TL 

Target  velocity  tables  in  launch 
coordinates 

zetgm(i) 

Kl 

IS1 

Print  flag,  IS1  f 1 indicates  print 
desired 

AAA 

a 

s 

Speed  of  sound 

111 

ha 

asL 

Specific  target  altitude  in  earth 
fixed  coordinates  + altitude  of 
missile  launch  site  above  sea  level 
(FT) 

PHO 

P 

Air  density  (SLUGS/FT3) 

QM 

Sn 

Mach  number  of  target 

WW 

vT 

Velocity  of  target  (Computed  in 

MAIN) 

CLA 

CL 

°T 

Aerodynamic  lift  coefficient  of 
target  due  to  angle  of  attack 

CLAA 

Table  containing  CLA  as  a function 
of  mach  number 

WT 

wT 

Weight  of  target 

ST 

ST 

Reference  surface  area  of  target 

XEU 

YEU 

Y^(°) 

Initial  acceleration  of  target  in 
earth  fixed  coordinates 

ZEU 

zTe(°) 

AT 

a 

Initial  target  angle  of  attack 

STL 


sin(e^) 


PROGRAM  VARIABLE 


SYMBOL 


DESCRIPTION 


Cl 

cos  (a) 

S9 

•'xfgCO)  + Y^CO) 

SI 

sin  (a) 

S2 

cos  (6^)  cos  (\pL) 

S4 

cos  (eL)  sin  (^) 

RLB 

l/B 

Apparent  ratio  of  plume  length  to 
breadth 

RLBK 

L/B 

Actual  ratio  of  plume  length  to 
breadth 

XLOS 

L 

X 

YLOS 

L 

Initial  components  of  T in  earth 

y 

fixed  coordinates 

ZLOS 

L 

z 

**' 

R 

R. 

1 

Initial  range  of  target  (FT) 

BET 

0 

cos  1 (y-qT!* 1 

L 

L 

Components  of  Horizon  Reference 

X 

X 

vector  in  earth  fixed  coordinates 

L 

L 

y 

y 

CLX 

% 

Components  of  target  center  line 

CLY 

% 

vector  in  earth  fixed  coordinates 

az 

TRP 

Trp 

Initial  plume  rotation  angle 

PROGRAM  VARIABLE 


SYMBOL 


DESCRIPTION 


THE 

e 

HC7FT 

Hc7 

Commanded  irradiance  for  plume 

P1IRSS 

Pi 

Vertex  angle  of  plume  transparency 

(DEG) 

P2IRSS 

P2 

Length  of  plume  transparency  (IN) 

AT7 

At7 

Area  of  plume  transparency 

RN 

A7 

Iris  ratio  number  7 

f 

Focal  length  of  projection  lens 

l 

Apparent  length  of  plume 

APPENDIX  A j 

This  Appendix  contains  a proposed,  structured  version  of  the  MICOM 
Hybrid  STINGER  program.  Usage  would  require  the  following  changes: 

t 

1.  Under  "Manual  Selection  of  Launch  Conditions"  where  the  variables 
KA,  ITT,  TAU,  IRSS  are  read  from  a card  - the  program  should  be 
altered  to  read  these  variables  from  the  terminal. 


2.  In  order  to  integrate  the  new  FLIGHT  and  REALT  programs  into  the 
structured  STINGER  package,  a few  changes  in  the  common  blocks  and 
variables  in  FLIGHT  and  REALT  will  be  necessary. 


3.  In  the  structured  version  of  STINGER  only  the  programs  MAIN  and 
INIT  have  been  essentially  altered.  Some  cleaning  up  and  alteration  may 
be  necessary  in  the  subroutines  which  perform  the  post  real-time  calcu- 
lations. 


4.  There  are  now  two  interpolation  routines  - INTERP  and  LINTRP. 
LINTRP  is  much  faster  and  INTERP  should  probably  be  removed  unless 
higher  order  interpolation  than  linear  is  desired. 


A.O 


J 


INPUTS  TO  STRUCTURED  STINGER  PROGRAM 


IS1 


NRUNS 


PRINT  FLAG 

IS1  = 0 for  limited  printout 
IS1  = 1 for  full  printout 
format  15 

DESIRED  NUMBER  OF  RUNS  PER  FLIGHT  CONDITION 
1 NRUNS  20 
format  15 


NPTS,  DX,  XTERM 

NPTS  = NUMBER  OF  DATA  POINTS  TO  BE  COLLECTED 
1 NPTS  50 

DX  = DELTA  RANGE  BETWEEN  COLLECTED  DATA  POINTS 

XTERM  = X SEPARATION  AT  WHICH  REAL  TIME  IS  TERMINATED 
format  110,  2F10,  0 


NL 


NUMBER  OF  SETS  OF  LAUNCH  CONDITIONS 
1 NL  61 
format  15 


(ZTEO(I),  ASC(I),  RR(I),  I = 1,  NL) 

THESE  ARE  THE  NL  SETS  OF  LAUNCH  CONDITIONS 

ZTEO(I)  IS  THE  I-TH  INITIAL  ALTITUDE 

ASC(I)  IS  THE  I-TH  CROSSING  ANGLE 

RR(I)  IS  THE  I-TH  INITIAL  RiN  GE 
format  3F10.  2 

III  LAUNCH  SELECTION  FLAG 

III  = 0 for  manual  selection 
III  = 1 for  automatic  selection 
format  15 

NOTE:  FROM  THIS  POINT  ON,  THE  INPUTS  TO  BE  SUPPLIED  DEPEND 

ON  THE  VALUE  OF  III  CHOSEN  ABOVE, 


THE  FOLLOWING  INPUTS  ARE  READ  ONLY  IF  III  = 0 (MANUAL  SELECTION) 
OR  AFTER  A FULL  SET  OF  AUTOMATIC  RUNS  IS  COMPLETED. 


A.l 


I\.n,  111, 


inu,  1 1\^£> 


KA  = NUMBER  OF  LAUNCH  SET  TO  BE  USED 
ITT  = NUMBER  OF  TARGET  TRAJECTORY  TO  BE  USED 
TAU  = TIME  POINT  ON  THE  TARGET  TRAJECTORY  AT  WHICH 
THE  MISSILE  IS  LAUNCHED 
IRSS  IS  THE  IRSS  FLAG 
IRSS  = 0 NOT  AN  IRSS  RUN 
IRSS  = 1,  IRSS  RUN 
format  215,  F5.  1,  15 

AFTER  A MANUAL  RUN  IS  COMPLETED,  THE  PROGRAM  WILL  CYCLE 
AND  READ  NEW  VALUES  OF  KA,  ITT,  TAU,  IRSS 

THERE  ARE  TWO  SPECIAL  CONVENTIONS  ON  KA. 

KA  negative  - restart  program  at  very  beginning 
KA  > 61  - terminate  program 

i 

THE  FOLLOWING  INPUTS  ARE  READ  ONLY  IF  III  = 1 (AUTOMATIC  SELECTION) 

NSAM  QUALITY  CONTROL  TEST  CASE  SAMPLE  SIZE 

NSAM-11 
format  15 

ICON  CONTINUATION  FLAG 

ICON  = 1 continue  after  quality  control  failure 
ICON  = 2 rerun  previous  runs  after  quality  control  failure 
format  15 

INL  INITIAL  LAUNCH  SET  NUMBER 

1 < INL  INL 
format  15 

(CLMS(l),  CLSl(I),  CLS2(I),  I = 1,  10) 

DATA  FOR  COMPUTING  CONTROL  LIMITS 
format  3 F6.  3 

Ql.  Q2,  Q3,  Q4,  Q5,  Q6 


DATA  FOR  COMPUTING  CONTROL  LIMITS 
format  6F10.  4 


VARIABLE  REFERENCE  LIST 


VARIABLE 

COMMON 

BLOCK 

USE 

1 

AAA 

AARUNS 

Speed  of  sound 

I 

AASUM 

A BRUNS 
ABSUM 

1 Averages  related  to  kill  probabilities 

ADSF(IO) 

COMM 

Array.of  analog-to-digital  scale  factors 

ARB(20) 

Array  for  storing  miss  distances  relative 
to  tailpipe 

ARBT(20) 

Array  for  storing  miss  distances  relative  tc 

ARC(20) 

See  ARB 

ARCT(20) 

See  ARBT 

ARM(20) 

See  ARB 

ARMT(20) 

See  ARBT 

AS 

EXTRA 

Target  crossing  angle 

ASC(61) 

Table  of  target  crossing  angles 

ASEO 

Target  dive  angle  (degrees) 

AT 

GTARG 

Initial  target  angle  of  attack 

AXB(20) 

See  ARB 

AXBT(20) 

See  ARBT 

AXC(20) 

See  ARB 

AXCT(20) 

See  ARBT 

AXM(20) 

See  ARB 

AXMT(20) 

See  ARBT 

A YB(20) 

See  ARB 

AYBT(20) 

See  ARBT 

AYC(20) 

See  ARB 

A YCT(20) 

See  ARBT 

A YM(20) 

See  ARB 

A YMT(20) 

See  ARBT 

AZB(20) 

See  ARB 

AZBT(20) 

See  ARBT 

AZC(20) 

See  ARB 

AZCT(20) 

See  ARBT 

A ZM(20) 

See  ARB 

AZMT(20) 

See  ARBT 

A 1 

Target  angle  of  attack  (RADIANS) 

A1P1 

COMP 

1 

A1P2 

COMP 

L 

A 1P3 

COMP 

I Control  limits 

A1P4 

COMP 

) 

A. 3 
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VARIABLE  REFERENCE  LIST  (Con't)  | 

1 

VARIABLE 

COMMON 

USE 

BLOCK 

i I 

A 1P5 

COMP 

1 

A1P6 

COMP 

1 Control  limits 

A2P1 

COMP 

A2P2 

COMP 

► 

A2P3 

COMP 

A2P4 

COMP 

A2P5 

COMP 

A2P6 

COMP 

BDHX(22) 

PK 

I 

BDHY(22) 

PK 

f Tables  related  to  lethality 

BDHZ(22) 

PK 

J 

CA 

Cosine  of  target  angle  of  attack 

CLA 

Target  lift  coefficient 

CLAA(IO) 

RTCOM 

Target  lift  coefficient  table  (function  of  mach) 

CLMS(IO) 

COMP 

1 

CLSl(lO) 

COMP 

l Data  for  calculating  control  limits 

CLS2(10) 

COMP 

! 

CLX 

) Components  of  initial  target  center  line  vector 

CLY 

V 

CLZ 

J 

COSE 

Cosine  of  the  angle  between  the  line  of  sight 

and  the  target  center  line 

COSEI 

RTCOM 

Initial  value  of  cose 

COSMU1 

MU 

! Terms  of  a transformation  matrix  (see 

COSMU2 

MU 

f subroutine  PRIME) 

COSMU3 

MU 

J 

COSSCO 

ANG 

Cosine  of  PSILR 

cosso 

ANG 

Cosine  of  THETLR 

CPL 

RTCOM 

= COSSCO 

CTL 

RTCOM 

= COSSO 

CX 

Components  of  initial  horizon  reference  vector 

CY 

Cl 

Cosine  of  AT 

DAOUT1 

-DA  Cl 

DAOUT2 

*DAC1 

DAOUT3 

-DA  Cl 

Variables  which  are  converted  from  digital 

DA OUT 4 

-DA  Cl 

to  analog 

DA OUT 5 

-DAC1 

DAOUT6 

-DAC1 

. 

A. 
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VARIABLE  REFERENCE  LIST  (Con’t.  ) 


VARIABLE 

COMMON 

BLOCK 

USE 

i 

DA OUT 7 

*DAC1 

DAOUT8 

*DAC1 

Variables  which  are  converted  from  digital 

DAOUT9 

*DAC1 

‘ to  analog 

DA  OUT  10 

*DAC1 

DA  OUT  11 

-DA  Cl 

DASF(ll) 

COMM 

Array  of  digital  to  analog  scale  factors 

DELTAR(30) 

MISL 

Target  range  table  (function  of  time) 

DT 

Current  time  in  real  time 

DTF 

RTC'OM 

Value  of  DT  at  termination  of  real  time 

DTR 

Conversion  factor-degrees  to  RADIANS 

DX 

EXTRA 

Delta  range  between  collected  data  points 

DX 

X separation  between  target  and  missile 

DXG 

RTCOM 

Value  of  DX  at  termination  of  real  time 

DY 

Y separation  between  target  and  missile 

DYG 

RTCOM 

Value  of  DY  at  termination  of  real  time 

DZ 

Z separation  between  target  and  missile 

DZG 

RTCOM 

Value  of  DZ  at  termination  of  real  time 

ECOS(17) 

RTCOM 

Table  used  in  calculating  IRIS  ratio  7 

EDOT 

RTCOM 

GGG/G 

FCI77(9) 

RTCOM 

Table  used  in  calculating  IRIS  ratio  7 

FI7T(9) 

RTCOM 

Table  used  in  calculating  IRIS  ratio  7 

FXA 

GTARG 

FXB 

GTARG 

FXC 

GTARG 

FYA 

GTARG 

Terms  of  a transformation  matrix  (see 

FYB 

GTARG 

• subroutine  MISCOM) 

FYC 

GTARG 

FZA 

GTARG 

FZB 

GTARG 

FZC 

GTARG 

G 

RTCOM 

Estimated  time  at  1000  feet  to  go 

GGG 

RTCOM 

Parameter  used  in  variable  analog  scale  fa 

HR7EE(17) 

RTCOM 

Table  used  in  calculating  IRIS  ratio  7 

IADIN1 

* A DC  1 

* 

IADIN2 

*ADC1 

1ADIN3 

-ADC1 

Values  resulting  from  analog  to  digital 

IADIN4 

-ADC1 

conversion 

LA  DIN  5 

*ADC1 

IADIN6 

*ADC1 

4 
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VARIABLE  REFERENCE  LIST  (Con't. ) 


VARIABLE 

COMMON 

BLOCK 

IADIN7 

*ADC1 

IADIN8 

*ADC1 

IADIN9 

-ADC1 

IADIN1G 

* A DC  1 

ICN 

MISD 

ICON 

IERR 

IFT 

COMP 

III 

MANEUV 

IIN 

*IDIS2 

INFLRT 

COMM 

INL 

lOUT 

*ODIS2 

IPASS 

IPLANE 

COMP 

IPTS 

RTCOM 

IQC1 

IRSS 

IS1 

PRNT 

ITAU(60) 

MANEUV 

ITRA  J(60) 

MANEUV 

ITT 

MANEUV 

188 

199 

COMP 

KA 

USE 


Values  resulting  from  analog  to  digital 
conversion 


Number  of  the  current  real  time  run 
(1  ICN  NR) 

Flag  to  indicate  continuation  after  quality 
control  failure 

Error  flag  for  equipment  reservation 
Flag  for  quality  control  failure 
Flag  to  indicate  manual  or  automatic  launch 
selection 

Input  discrete  word 

Flag  to  indicate  first  pass  through  real  time  code 
Launch  set  number  (automatic  selection) 

Output  discrete  word 

Flag  to  indicate  first  pass  through  main  program 
Parameter  used  in  quality  control  test 
Number  of  data  points  to  be  collected  before 
intercept 

Flag  to  indicate  first  launch  set  (automatic  sel. ) 

Flag  to  indicate  IRSS  run 

Flag  to  indicate  full  or  partial  printout 

Table  of  launch  times  (see  TAU) 

Table  of  target  trajectory  numbers 
Number  of  current  target  trajectory 
Flag  to  indicate  last  launch  set  (automatic 
selection) 

Flag  to  indicate  test  case  is  due 
Number  of  current  launch  conditions 


Table  related  to  lethality 

Flag  to  indicate  target  maneuver  in  progress 
Counters  for  the  two  lodps  in  the  real  time 
code 

Maximum  execution  time  of  real  time  loop 
Number  of  entries  in  range  table  (see  XXS) 
Error  flag  for  subroutine  INTERP 
Number  of  sets  of  launch  conditions 


■- ^ - 


VARIABLE  REFERENCE  LIST  (Con't.  ) 


VARIABLE 

COMMOl 

BLOCK 

NPATH 

COMM 

NPTS 

EXTRA 

NPX 

NR 

NRUNS 

NSAM 

COMP 

NT 

RTCOM 

NX 

PHO 

PI 

PI02 

PPX(50) 

RTCOM 

PFY(50) 

RTCOM 

PPZ(50) 

RTCOM 

PSIL 

PSILR 

QM 

QMM(IO) 

RTCOM 

Ql 

COMP 

Q2 

COMP 

Q3 

COMP 

Q4 

COMP 

Q5 

COMP 

Q6 

COMP 

RI 

RTCOM 

RLB 

RLBI 

RTCOM 

RLBK 

RTCOM 

RM(3) 

MISD 

RN 

RNI 

RTCOM 

USE 


Flag  to  indicate  which  loop  to  execute  in 
real  time  code 

Desired  number  of  data  points 
Number  of  data  points  used  in  interpolation 
(NPX  = 2 is  linear  interpolation) 

Number  of  flights  to  be  performed  in  current 
run 

Desired  number  of  runs  per  flight  condition 
Quality  control  test  case  sample  size 
(2  NSAM  11) 

Number  of  points  in  target  trajectory  tables 
Argument  for  subroutine  interp-number  of 
points  in  tables 
Air  density 
3. 14159265 
PI/ 2. 

Arrays  for  storing  DX,  DY,  and  DZ 

respectively  at  points  corresponding 
to  range  table  entries  (see  XXS) 

PSI  launch-one  of  the  two  euler  angles  in  the 
transformation  from  earth  to  launch 
coordinate  system  (degrees) 

PSIL  converted  to  RADIANS 
Target  mach  number 

Table  of  mach  numbers  (goes  with  CLAA(IO)) 


Data  for  calculating  control  limits 


Initial  range-missile  to  target 

Apparent  plume  length  to  breadth  ratio 

Initial  value  of  RLB 

Actual  plume  length  to  breadth  ratio 

MISS  distance  array 

IRIS  ratio  7 

Initial  value  of  RN 
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VARIABLE  REFERENCE  LIST  (Con't.  ) 


r 


* 


VARIABLE 

COMMON 

BLOCK 

USE 

l 

RR(61) 

Table  of  target  initial  ranges 

RTD 

Conversion  factor-radians  to  degrees 

SA 

Sine  of  target  angle  of  attack 

SEO 

Target  dive  angle  (RADIANS) 

SINMU1 

MU 

| Terms  of  a transformation  matrix  (see 

SINMU2 

MU 

r subroutine  PRIME) 

SINMU3 

MU 

J 

SINSCO 

ANG 

Sine  of  PSILR 

SINSO 

ANG 

Sine  of  THETLR 

SLO 

Sense  line  zero 

SL5 

Sense  line  five 

SL6 

Sense  line  six 

SL15 

Sense  line  fifteen 

SPHG 

PK 

Variables  related  to  kill  probabilities 

SPHS 

PK 

SPK 

PK 

SPL 

RTCOM 

= SINSCO 

SPM 

PK 

See  SPHG 

SPO 

Plume  rotation  angle 

SPOI 

RTCOM 

Initial  value  of  SPO 

STL 

RTCOM 

= SINSO 

SI 

Sine  of  AT 

S1P1 

COMP 

S1P2 

COMP 

Control  limits 

S1P3 

COMP 

S1P4 

COMP 

S1P5 

COMP 

S1P6 

COMP 

S2 

RTCOM 

Term  of  earth- to-launch  transformation  matrix 

S2P1 

COMP 

’ 

S2P2 

COMP 

S2P3 

COMP 

» Control  limits 

S2P4 

COMP 

S2P5 

COMP 

S2P6 

COMP 

S3 

RTCOM 

\ 

S4 

RTCOM 

V See  S2 

S5 

RTCOM 

J 

8 


VARIABLE  REFERENCE  LIST  (Con't. ) 


r 


VARIABLE 

COMMON 

BLOCK 

TAMA(30) 

MISL 

TAU 

THETAL 

RTCOM 

THETLR 

TIME(50) 

RTCOM 

TM(30) 

TMA(30) 

RTCOM 

TRP 

VM(30) 

MISL 

VMX(50) 

RTCOM 

VMY(50) 

RTCOM 

VMZ(50) 

RTCOM 

VT 

GTARG 

VTI 

VVVV 

EXTRA 

W1 

COMP 

W2 

COMP 

W3 

COMP 

W4 

COMP 

W5 

COMP 

W6 

COMP 

H 

X 

XC 

XCOMP 

XD 

COMB 

XDDE 

COMB 

XDEF 

COMB 

XDEM(30,  45) 

MANEUV 

XDH(22) 

PK 

XDM(30) 

RTCOM 

XDO 

RTCOM 

' 

USE 


Table  of  times-goes  with  DELTAR  (30)  and 
VM  (30) 

Time  point  on  target  trajectory  at  which 
missile  is  launched 
THETA  launch  - see  PSIL  (degrees) 

THETAL  converted  to  RADIANS 
Array  for  storing  DT  at  points  corresponding 
to  range  table  entries  (See  XXS) 

Table  of  times-goes  with  XDEM,  YDEM,  ZDEM 
Table  of  times-goes  with  XDM,  YDM,  ZDM 
= SPO 

Target  velocity  table  (function  of  time) 

Arrays  for  storing  XDOT,  YDOT  and  ZDOT 
Respectively  at  points  corresponding 
To  range  table  entries  (see  XXS) 

Target  velocity  (in  subroutine  MISCOM) 

Target  velocity  in  real  time 
Initial  target  velocity 

Data  for  calculating  control  limits 


Initial  target  X corrdinate  in  earth  system 

Target  velocity  X component- earth  system 

Target  velocity  X component-launch  system 

Used  in  earth-to-launch  transformation 
(subroutine  GUID) 

Target  acceleration  over  data  collection 
range-X  component 

Target  velocity  at  first  range  table  entry-X 
component 

Array  of  45  target  velocity  tables  in  earth 
system-X  components 

Table  related  to  lethality 

Current  target  velocity  table  in  earth  system-X 
component 

Value  of  XDOT  at  termination  of  real  time 


VARIABLE  REFERENCE  LIST  (Con't.  ) 


VARIABLE 

COMMON 

BLOCK 

USE 

i 

XDOT 

Target-missile  relative  velocity-X  component 

XDTG 

COMB 

See  XD 

XDTGMS(30) 

RTCOM 

Current  target  velocity  table  in  launch 
system-X  component 

XDTGO 

RTCOM 

Target  initial  velocity  in  launch  system-X 
component 

XE(30) 

RTCOM 

Target  acceleration  table-X  component 

XLOS 

Missile  to  target  line  of  sight-X  component 

XM(3) 

MISD 

MISS  distance  array 

XMAN14,  50) 

RTCOM 

Array  for  storing  DT,  XCOMP,  YCOMP, 

ZCOMP  at  range  table  entries 

XMISS(7) 

RTCOM 

Array  for  storing  DX,  DY,  DZ,  XDOT,  YDOT, 

ZDOT  and  XMISSC  if  missile  misses 
target 

XMISSC 

Variable  which  indicates  a miss 

XTERM 

EXTRA 

X separation  between  missile  and  target  at 
which  real  time  is  terminated 

XXS(50) 

RTCOM 

Range  table  for  data  collection-NADJ  points 
separated  by  DX  and  ending  with  XTERM 

XXX 

X separation  between  target  and  missile- 
comes  from  an  analog  value  which  is 
scaled  ten  times  smaller  than  DX 

Y 

Initial  target  Y coordinate  in  earth  system 

YA 

COMD 

Related  to  MISS  distance  1 

YC 

Target  velocity  Y component- earth  system 

YCOMP 

Target  velocity  Y component-launch  system 

YD 

COMB 

See  XD 

YDDE 

COMB 

See  XDDE-Y  component 

YDEF 

COMB 

See  XDEF-Y  component 

YDEM(30,  45) 

MANEUV 

See  XDEM-Y  components 

YDH(22) 

PK 

See  XDH 

YDM(30) 

RTCOM 

See  XDM-Y  component 

YDO 

RTCOM 

Value  of  YDOT  at  termination  of  real  time 

YDOT 

See  XDOT-Y  component 

YDTG 

COMB 

See  XD 

YDTGMS(30) 

RTCOM 

See  XDTGMS-Y  component 

YDTGO 

RTCOM 

See  XDTGO- Y component  ! 

YE(30) 

RTCOM 

See  XE-Y  component  J 

VARIABLE 


YLOS 

YM(3) 

YYY 

Z 

ZA 

ZABOVE 

ZALT 

ZC 

ZCOMP 

ZD 

ZDDE 

ZDEF 

ZDEM(30,  45) 

ZDH(22) 

ZDM(30) 

ZDO 

ZDOT 

ZDTG 

ZDTGMS(30) 

ZDTGO 


ZTEO(61) 

ZZZ 


VARIABLE  REFERENCE  LIST  (Con't.  ) 


COMMON 

BLOCK 


COMD 


RTCOM 


COMB 

COMB 

COMB 

MANEUV 

PK 

RTCOM 

RTCOM 

COMB 

RTCOM 

RTCOM 

RTCOM 


See  XLOS-Y  component 
MISS  distance  array 
See  XXX-Y  component 

Initial  target  Z coordinate  in  earth  system 
See  YA 

Altitude  of  launch  site  (negative) 

Target  altitude  (negative) 

Target  velocity  Z component-earth  system 
Target  velocity  Z component-launch  system 
See  XD 

See  XDDE-Z  component 
See  XDEF-Z  component 
See  XDEM-Z  component 
See  XDH 

See  XDM-Z  component 

Value  of  ZDOT  at  termination  of  real  time 
See  XDOT-Z  component 
See  XD 

See  XDTGMS-Z  component 

See  XDTGO-Z  component 

See  XE-Z  component 

See  XLOS-Z  component 

MISS  distance  array 

Table  of  target  initial  altitudes 

See  XXX- Z component 


COMMON  BLOCK  CROSS-REFERENCE  TABLE 


ANG 

COMB 

COMD 

COMM 

COMP 

EXTRA 

GTARG 

MANEUV 

MISD 

MISL 

MU 

PK 

PRNT 

2 
o 
u 
t-> 
g • 

MAIN 

X 

1 

X 

X 

X 

X 

X 

1 

X 

X 

X 

BLOCK  DATA 

■ 

■ 

X 

X 

X 

X 

FLIGHT 

■ 

X 

■ 

■ 

X 

GUID 

X 

X 

■ 

■ 

HWFLOT 

■ 

■ 

■ 

■ 

INIT 

X 

* 

X 

INTERP 

■ 

■ 

K SCALE 

■ 

■ 

X 

X 

X 

LETH 

■ 

X 

X 

X 

X 

X 

X 

X 

LINTRP 

■ 

■ 

■ 

■ 

■ 

MEAN 

■ 

■ 

X 

■ 

■ 

■ 

■ 

■ 

MISCOM 

X 

X 

■ 

X 

X 

X 

X 

X 

X 

PRIME 

X 

■ 

■ 

■ 

■ 

X 

X 

X 

X 

X 

QCLIM 

■ 

■ 

X 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

REALT 

■ 

X 

■ 

■ 

■ 

■ 

■ 

■ 

■ 
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cosscc 

cosso 

SIN  SCO 
SINSO 


COMMON/COMB/ 


XD 

XDDE 

XDEF 

XDTG 

YD 

YDDE 

YDEF 

YDTG 

ZD 

ZDDE 

ZDEF 

ZDTG 
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COMMON/COMD/ 


YA 

ZA 


COMMON/COMM/ 


ADSF  (10) 

DASF  (11) 

INFLRT 

LOOP1 

LOOP2 

LTMAX 

npath 


A. 16 


A1P1 

A1P2 

A1P3 

A1P4 

A1P5 

A1P6 

A2P1 

A2P2 

A2P3 

A2P4 

A2P5 

A2P6 

CLMS  (10) 

CLSl(lO) 

CLS2(10) 

IFT 

IPLANE 

199 

NSAM 

Ql 

Q2 

Q3 

Q4 

Q5 

Q6 


SI  P2 


COM  MON /EXTRA/ 


AS 

DX 

NPTS 

VVVV 

XTERM 


COMMON/  MA  NEUV/ 


B 

hi 

ITAU(60) 
ITRA  J(60) 
ITT 
NADJ 

XDEM(30,  45) 
YDEM(3 0,  45) 
ZDEM(30,45) 


-i 

i 


r 


ICN 

RM(3) 

XM(3) 

YM(3) 

ZM(3) 


COMMON/ MISD/ 


COMMON/ MISL/ 


DELTAR(30) 

TAMA(30) 

VM(30) 


A, 22 


- 


- 


COMMON/MU/ 


COSMU1 

COSMU2 

COSMU3 

SINMU1 

SINMU2 

SINMU3 


m 


COMMON/PK/ 

BDHX(22) 

BDH7(22) 

BDHZ(22) 

KILU22) 

SPHG 

SPHS 

SPK 

SPM 

XDH(22) 

YDH(22) 

ZDH(22) 


A.  24 


COMMON/RTCOM/ 


CLAA(IO) 

COSEI 

CPL 

CTL 

DTF 

DXG 

DYG 

DZG 

ECOS(17) 

EDOT 

FCI77(9) 

FI7T(9) 

G 

GGG 

HR7EE(17) 

IPTS 

LEVEL 

NT 

PPX(50) 

PPY(50) 

PPZ(50) 

RI 


RLBI 

RLBK 

RNI 

SPL 

SFOI 

STL 

52 

53 

54 

55 

THETA L 

TIME(50) 

TMA(30) 

VMX(50) 

VMY(50) 

VMZ(50) 

XDM(30) 

XDO 

XDTGMS(30) 
XDTGO 
XE(30) 
XMAN(4,  50) 
XMISS(7) 
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r 

h 


COMMON/RTCOM/  (Con't) 

XXS(50) 

I 

YDM(30) 

YDO 

YDTGMS(30) 

n 

YDTGO 

YE(30) 

ZALT 

ZDM(30) 

ZDO 

ZDTGMS(30) 

ZDTGO 

ZE(30) 
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EXPLANATION  OF  FLOW  DIAGRAMS 


The  flow  diagrams  included  here  are  hierarchical  in  nature  with  a 
top-down  structure.  In  the  first-level  flow  diagram,  the  program  is  broken 
down  into  its  major  modules,  numbered  I. , II.  , III. , etc. , and  the  relation- 
ships among  these  modules  are  shown. 

Each  module  in  the  first-level  diagram  is  analyzed  in  more  detail 
in  a second-level  flow  diagram.  For  example,  module  III.  in  the  first-level 
diagram  is  broken  down  into  submodules  numbered  III.  A.  , III.  B.  , III.  G.  , 
etc. 


Some  of  the  second-level  modules  are  broken  down  further  in  third- 
level  diagrams.  For  example,  module  III.  K.  in  the  second-level  diagram 
devoted  to  module  III.  is  broken  down  in  a third-level  diagram  into  blocks 
numbered  III.  K.  1. , III.  K.  2. , etc. 

The  symbol  + appearing  before  a block  number  indicates  that  there 
is  a lower-level  diagram  devoted  to  that  block. 


MAIN  PROGRAM 

STINGER  HYBRID  SIMULATION 

FIRST -LEVEL  FLOW  DIAGRAM 


1 


: 

ii 


MAIN  PROGRAM 

STINGER  HYBRID  SIMULATION 

SECOND-LEVEL  FLOW  DIAGRAMS 


A, 31 


j 


V. 

SELECTION  OF  TEST 
CASE  LAUNCH  CONDITIONS 


A.  37 


r 


r 


MAIN  PROGRAM 

STINGER  HYBRID  SIMULATION 


THIRD  LEVEL  FLOW  DIAGRAMS 


ADJUST  TARGET  VELOCITY 
TABLE  FOR  LAUNCH  TIME 


CALCULATE  ANGLE  BETWEEN 
POSITIVE  X-EARTH  AXIS  AND 
TARGET  INITIAL  VELOCITY  VECTOR 


VI. D. 3 


VI. D. 4 


VI. D. 5 


VI. D. 6 


VI. D. 7 


VII. B. 

CLEAR  STORAGE  AREAS 


VII. E. 

CHECK  TIME  CRITICAL 
DATA  COLLECTION 


VII. I. 

CALCULATE,  SAVE  AND  PRINT  MISS 
DISTANCE  DATA  REFERENCED  TO  TAILPIPE 


VII. I. 1. 


VII. I. 2. 


VII. I. 3. 


VII. I. 4. 


VII. I. 5. 


VII. K. 

CALCULATE  AND  SAVE  MISS  DISTANCE 
DATA  REFERENCED  TO  CG 


VII. K.l. 


<=> 


AXMT(ICN)  = XM (1) 
AYMT(ICN)  = YM (1) +19 . 5*YA 
AZMT(ICN)  = ZM (1) +19 . 5*ZA 


VII. K. 4. 


VII. K. 5. 


ARBT(ICN)  = SQRT  „((XM  (2)  -19 . 5)  **2+YM  (2)  **2 
+ZM (2) **2) 


VII. K. 6. 


ARCT(ICN)  = SQRT  ;((XM  (3)  -19 . 5)  * *2+YM  (3) **2 

+ZM (3) **2) 


,5 


VIII. D. 

DECIDE  WHERE  TO  TRANSFER  CONTROL 
NEXT  BASED  ON  MANUAL  OR  AUTOMATIC 
LAUNCH  SELECTION,  TEST  CASE  RUN,  QUALITY 
CONTROL  FAILURE,  ETC. 


ENTER 


VIII. D.l. 


QUALITY^. 
CONTROL  FAILURE^* 
I FT  = 1? 


VIII. D. 2. 


CASE  RUN? 
199  = 1?^ 


VIII. D. 3. 


X^AUTOMAT  ICS. 
WINCH  SELECTION^ 
^ III  = 1? 


VIII. D. 4 


VIII. D. 5. 


199  = 0 


■ "^FiRST^, 
NO  <^TEST  CASE? 
^\IQC1  = 1?, 


VIII. D. 7.  JL 

N°XtEST  CASE?X.  - YES 
^'V'nI88  = 


VIII. D. 8.  I88  = I 


, VIII. D. 9.1 


next  page 


VIII. D. 6.  llQCl  = 0 


VIII. D.  (Continued) 
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UNCLASSIFIED 

3 If. 
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V. 

CALCULATE  RANGE  TABLE 
FOR  DATA  COLLECTION 


: 


A.  68 


CALCULATE  INITIAL  COSE 


VII. C. 


XI. 

CALCULATE  INITIAL  TARGET 
CENTER  LINE  VECTOR 


calculate  initial  plume 

ROTATION  ANGLE 


XII. A. 


XII. B. 


XII. C. 


XII. D. 


XII. E. 


XII. F. 


XII. G. 


5 


CALCULATE  INITIAL 
IRIS  RATIO  7 


ENTER 


XIII. A. 


NX  = 17 
NPX  = 2 


XIII.B.  I CALL  INTERP (COSEI , ECOS,  HR7EE,  NX,  NPX,  HR7E,  NERR) 


XIII. C.  I HC7FT  - 1209 . 675*HR7E*EXP (-2 . 341*ALOG (RI/22965 . 831) ) 


XIII. D. 


XIII. E. 


XIII. F. 


XIII. G. 


P1IRSS  = 2 . *ATAN ( . 5/RLBI ) 


P2IRSS  = 22 . 5*RLBI/RI 


AT7  = ( (P2IRSS**2) *SIN (P1IRSS) ) / 
(2.* (COS (P1IRSS/2 .) )**2) 


PJTU1  = .1*EXP (1 . 003258*ALOG (AT7) ) 


XIII. H. 


XIII. I. 


FI 7 = HC7FT/PJTU1 


NX  = 9 
NPX  = 2 


XIII. J. 


CALL  INTERP (FI7,  FI7T,  FCI77,  NX,  NPX,  RNI , NERR) 


RETURN 
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THUNK  LINE  PATCHING  TS 


r 


* 


c 

c 

c 

c 


c 

c 

c 


v> 1 


PROGRAM  TO  CHECK  DADIOS  DISCRETES*  THIS  IS  ACCOMPLISHED  BY  PATCHING 
INPUT  AND  OUTPUT  DISCRETE  WORDS  TOGETHER  AT  THF  AD/A 1 LOGIC  PANNEL 


THE  OUTPUT  DISCRETE  COUNTS  UP 
BEGINS  OVER  WITH  0. 


NAAGS»T10»CM70000»PT  • 

ACCT  (PN=PARKERL  »PBC=7487G1 3R00  *CC=7300 ♦ 0P=A3*  JN= J002) 
fTNH.  ■ . ' 

SEGLOAD.  V ' • . ■ . ; • ■ 

EGO.  ' ' " ' \ •' 

0000000000000000000000  v . 

PROGRAM  TRDISIO (OUTPUT. HFILE*TAPE6=0UTPUT) 

C 

c 

c , 
c 
c 
c 
c 
c 

C-  ' 

c 

c . 

C 

C ; 

c 

C > 

c 

C 
C 

c 

C 
C 


CYCLICLY  FROM  0 TO  2*«16-1  AND  THEN 


PROGRAM 
I ERR 

I ST  AT 


IIDIS 

IODIS  : 

MAXLINE 

MAXPAGE 

IOUT 

LEN 

,PIN(I) 
POUT ( I ) 


VARIABLES 


ERROR  CODE  FOR  RESERVATION  . : 

0=NOERROR,  GT.O=RESERVATTON  ERROR 
REAL  TIME  MODE  ^ / ' ■ / 

0=IN  REAL  TIME»  (ISTAT.GT.O  NOT  ;IN  REAL  TIME 
CDC-6600  SENSE  LINE  DISCRETE <16_BIT>  IIDIS=IBACK 
CDC-6600  CONTROL  LINE  DISCRETE (16  BIT)  JODIS=IOUT 
MAXIMUM  LINES  PER -PAGE  ; / 

MAXIMUM  PAGES  OF  PRINTOUT  , Y-  - ' ' 

A COUNTER  WHICH  PRODUCES  OUTPUT  FOR  TESTING 
NUMBER  OF  TRUNKING  STATION  CONNECTOR  PAIRS  USED 
I-TH  'TRUNCKING  STATION  CONNECTOR  INPUT  - ? 

l-TH  TRUNCKING  STATION  CONNECTOR  OUTPUT  - 


DADIOS  PATCHING  REQUIREMENTS  (ONE  IIDIS  AND  ONE  IODIS) 


), 


, C 

..  c 

TRUNKING 

'•  \ : W .'V'  Y-  1 {■' 

FORTRAN  V •• 

••••  • • ■ y ■ ' y , ■ .■  • ’ !.'■ 

‘ AD/4 "PATCH  ' 

c 

V-50  TO  W-50 

’ FOR  ./IDISP/I .IIDIS 

TR0Q-TR07  AND  TR20-TR27 

c 

V-r 52  TO  W-50 

F OR ; /IDI S?/l * I IDI S 

YTR40-TR47  AND  TR60-TR67  * 

t C ■ 

V-50  TO  W-51  • 

FOR  /IDIS2/2. IIDIS  , 

TR00-TR07  AND  TR20-TR27 

. c ,*■ 

r 

V-52  TO  W-51  ' • 

FOR  . /IDIS?/P« I IDIS 
• ; 

, TR40-TR47  AND  TR60-TR67  ! 

V 

. c 

V-51  TO  W-60 

FOR  /0DIS2/1. IODIS 

TR10-TP17  AND  TR30-TR37 

c 

V-53  TO  W-60 

FOR  /0DIS2/1.I0DIS 

TR50-TR57  AND  TR70-TR77  . 

V-51  TO  W-61  ' FOR  /0DIS2/2. IODIS 
V-53  TO  W-61  . FOR  /0DIS2/2 » IODIS 


TRl 0-TR17 ; AND 
TR50-TR57  AND 


TR30-TR37 

TR70-TR77 


j ■ \ 


INTEGER  PIN"(8)  .POUT  (8) 

DIMENSION  IBITIN(16) ♦ IBITOUT (16) 
COMMON/INTCOM/IOUT .LINE .MAXL'INE .ISA VEIN (100) 
INTERRUPT (1=1 *R=50.T=500*P=1000) 
C0MM0N/*IDIS2/1*IIDIS  - , 

C0MM0N/*0DIS2/1. IODIS 

INITIALIZATION 


I0UT=0 

L 1 NE=0 

MAXLINE=60 

MAXPAGE=10 

IPAGE=0 

LEN=4 

PIN ( 1 ) =3RV50 
PIN  (2)  =3RV5.1 
PIN (3) =3RV52 


: > I -’ 


i 


J 


PIN(4)=3RV53 


CALL  RESERVE (IERR) 
WRITE (6*1000)IERR 
IF ( IERR«NE*0) STOP 


CHECK  PATCHING  STATUS 


CALL  PATSTAT (PlN*POUT*LEN) 
00  10  1=1 *LEN 

10  WRITE(6*8000)PIN(I> »POUT(I) 


REAL  TIME 


CALL  .SIMRUN  (ISTAT)  , - 

WRITE (6»2000) ISTAT 

CALL  REMARK  (17H  JOB  ,IM  REAL  TIME) 


25  CONTINUE 

CALL  BHOLD  • 

LINE=0 

WRITE (6*5000) 
IPAGE=IpAGE*l 
IF ( I PAGE •EQ.IOJGO  TO  20 
11=1 OUT-MAXLINE* 1 


DO  30  1=11* I OUT 
CALL  BITS(IMBITOUT) 

;LINE=LINE*1  " , - 

IO'JMY=ISAVEIN  (LINE) 

0 CALL  BITS ( IDUMY * IP IT IN) 

30  WRITE (6 *4000) I* IBITOUT* IDUMY* IBITIN 
L1NE=0  ’ 

CALL  SIMGO  )■ 

GO  TO  '25  ■ 

20  WRltE (6*3000) 

CALL  REMARK (15H  RETURN  TO  MAIN) 

STOP  , ' , , 


1000  F0RMAT(24H1RESERVATI0N  ERROR  C0DE=*0?0) 

2000  F0PMATU8H  PEAL  TIME  STATU5=*020) 

3000  FORMAT ( 1 HO ♦ ‘PROGRAM'  TERMINATED  NORMALLY*) 

4000  FORMAT (10X*I10*5X*16I1*SX*T10*5X*16I1) 

5000  FORMAT (1H1) 

*000  FORMAT (/*OTRUNK  LINE  PATCHING  IS  *,R3**  TO  **P3) 
END  ' • ' 

SUBROUTINE  SUB1 


REAL  TlMf  INTERPUPT  SUBROUTINE 


COMMON/INTCOM/IOUT,LINE*MAXLINE*ISAVEINdOO) 
COMMON/* ID1S2/1 * I IDIS 
COMMON/*ODIS2/1*IODIS  , 


on  o on  noon 


IOOIS=IOUT  * 

IN=IIDIS 

LINE=LINE*1 

ISA VEIN (LINE) =IN  ' 

IF  (I OUT •GE*2**16-1 ) I0UT=0 

IF ( I OUT  «LT ,2** 16-1 ) 1 0UT= I OUT ♦ 1 

IF(LINE.GE.MAXLINE)CALL  SImHOLD 

CALL  SIMTDLE 

END 

SUBROUTINE'  PITS  (I IN* IBIT) 

PPQGRAM  TO;CONVFPT  BASE  TEN  TO  BITS 


THE  LEAST  SIGNIFICANT  BIT  IS  STORED  IN  BIT(1) 


PROGRAM  varibles 
I IN 

IBIT  ( I ) / . 

MAXBIT - ■ 


A NUMBER  BASE  TEN  v V 

I-TH  BIT  OF  A BINARY  NUMBER  • 

MAXIMUM  BITS  CONVERTED  IN  INPUT  SENSE  LINE 


INTEGER  OLDNUM 
DIMENSION  IBIT (16) 

, MAXBIT=16 

DO  1 1=1 »MAXBIT 

1 IBIT ( I ) =0  , 

0LDNUM=I IN 

DO  2 1=1 *MAXBIT 

NEWNUM=0LDNUM/2 

IBIT (I ) =0LDNUM-2*NEWNUM 

OLDNUM=NEWNUM 

2 CONTINUE  , 

IHALF-MAXBIT/2 

DO  3 1=1 ♦IHALF  ; 
ITEMP=IBIT  ( I ) 

IBIT (I)=IBIT (MAXBIT*1-I) 

3 IBIT (MAXBIT+l-I ) =ITEMP 
RETURN 

END- 

0000000000000000000000  ) 

RTREE  TRDISIO(O) *SUB1 (1) 
, ‘ GLOBAL  INTCOM  v 

END 

oooooooooooooooooooood  • • ; 
0000000000000000000000 
oooooooooooooooooooood  «- 


ooo  o o o o o o r>  o o r>  o r>  o o r>  r>  o o o o o o r>  o o o o r>  o o o o o’  r>  o O o 


NFEGS.T25.CM70000.RT. 

ACCT(PN=PARKERL.PBC=74R7G13R00.CC=7300.0P=A3. JN=J0C?> 

FTNH.  .• 

SFGLOAD.  ‘ . , 

LGO.  - • - ' w ’ ...  v 

0000000000000000000006 

PROGRAM  TRALGIO (OUTPUT .HF IlE »TAPE6=0UTPUT ) 

•’  . ••  . * (’•  ’ • 

PROGRAM  TO  TEST  ANALOG  SIGNALS  BETWEEN  AD/A  AND  CDC/6600.  THIS  TASK 
IS  ACCOMPLISHED  BY  TURNING  THE  ANALOG  SIGNAL  AROUND  AT  THE  AO/4  AND 
COMPAREING  DIFFERENCE  UPON  RETURN  TO  THE  CDC/6600.  THE  PROGRAM  TEST 
FOR  ERRORS  GREATER^  THAN  FIVE  PERCENT. 

NOTE  ADCS  AND  DACS  ARE  IN  GROUPS  OF  16*  1-16,  17-32*  33-4B,  49-64. 
FLOATING  POINT  ANALOG  SIGNALS  APE  SCALED  GE  -l'.O  AND  LE  ♦ 1.6. 
INTEGER  ANALOG  SIGNALS  ARE  SCALED  GE  -32767  AND  LE  ♦32767(14  BIT). 

V . ‘ ■ . i:A  ...  • ' . ' • 

> PROGRAM  VARIABLES  ‘ ' : V t, 

IERR  ERROR  CODE  FOR  RESERVATION  . V \ • 

■ '•  0=N0ERR0R,  GT.O=RESERVATlON  ERROR 
ISTAT  •'  REAL  TIME  MODE  - : ; , 

0=IN  REAL  TIME*  ISTAT.GT.O  NOT  IN  REAL  TIME  " - 

OUT  -1'  V > THE  DAC  VARIABLE 
■ - BACK  THE  ADC  VARIABLE  \ 

V BIGDIFF  , MAXIMUM  ALLOWABLE  DIFFERENCE  TN  ADC  AND  DAC  . 

DIFF  ACTUAL  COMPUTED  DIFFERENCE  IN  ADC  AND  DAC 

LINE  NUMBER  OF  LINES  OF  PRINTOUT  IN  EXECUTION 

MAXLINE  MAXIMUM  NUMBER  OF  LINES  PRINTFD 

TIMELFT.  REMAINING  PCT  TIME  ' ’ ' ' 

STATUS  /-THE  HRTM  STATUS  WORD 

LOOP  NUMBER  OF  ATTEMPTS  TO  EQUATE  ADC  AND  DAC 

J A MASKING  WORD  USED  IN  COJUNCTION  WITH  THE 

, STAT'  FUNCTION  ’ - 

DADIOS  PATCHING  REQUIREMENTS  (AD/4  FIELD  3*  4TH  GROUP  ADC  AND  DAC) 

v - TRUNKING  • FORTRAN  AD/4  PATCH 

W-03  TO  V-06  FOR  /*ADC1 /49* ADC  ' TRIO-T R17  AND  -TR30-TR37 

W-l 3 TO  V-07  FOR  /*OAC 1/49. DAC  TR50-TR57  AND  TR70-TR77 

* • ' • . • * f 4 * • 

THE  FIRST  ADC  IN  EACH  GROUP  OF  16  CORRESPONDS  TO  - AD./ 4 TR-10. 

THE  FIRST  DAC  IN  EACH  GROUP  OF  16  CORRESPONDS  TO  AD/4  TP-50.  . 


COMMON/INTCOM/OUT*LOOP*BACK(16) .BIGDIFF *DIFF( 16) .STATUS* J 
INTERRUPT (1=1 .R=l 00, T= 10006) 

C0MM0N/*ADC1/  l.ADC(16>  ..  ' 

C0MM0N/*DAC1/  1, DAC (16)  ' f ' 

INITIALIZATION  . s \ 

* • . • \ i f ' * 

BIGDIFF*. 001  .*  .:  • 

OUT =-l  .5  ■ .'i  . 

LOOP=0 

L INE=0  ' . • 

MAXLINE=400  ' t ’ 


"WRITE (6* 1000 ) IERR 

IF(IERR*NE*0)ST0P  ' ' " . 

C ' ; . V ' • - ■ ' * • • 

C REAL  TIME  * y ./ 

* ,C  • • i ■ • ; if'  - i 

! ‘ CALL  SIMRUN(ISTAT) 

WRITE  (6*2000)  1ST  AT  v. 

CALL  REMARK (17H  JOB  IN  REAL  TIME) 

WRITE (6*5000) 

, 25  CONTINUE  v ' 

call  bhold  " 

5 LINE=LINE*1 

IF (LINE* GT» MAXLINE) GO  TO  50 

6 - -TIMELFT-STATUS.A.3777777B 

WP ITE (6 *4000 ) OUT *BACK*LOOP  ( 

}■  WRITE  (6 *6000)  DIFF 

. 0UT=0UT*0.01  , 

l >>  IF(0UT.GT.1.56)0UT=6.0 
L00P=0 

! CALL  SIMGO  / ' v 

■ I v GO  TO  25  • ;.>v  - \ . ■ - ... 

50  WRITE(6.3000) 

■ stop  ' ■£$  h:r-(  • ••  ’ ' t ' • - i. 

1000  FORMAT (24H1 RESERVATION  ERROR  CODE=*02O) 

••  2000  FORMAT  ( 18H  REAL  TIME  STATUS=*020), 

■ 3000  F0RMAT(1H0**PR0GRAM  TERMINATED  ON  MAX  LINF*) 
4000  F ORMAT ( IX* 17F7.4  * 13) 

5000  FORMAT(lHl)  . . ?r.( 

6000  FORMAT (8X*16F7*4)  . - - , 

END  ’ ' ' i . 

SUBROUTINE  SUB1  v ' . . 

C . • . • - 

C REAL  TIME  INTERRUPT  SUBROUTINE 


, COMMON/INTCOM/OUT, LOOP, BACK (16) «BIGDIFF*DIFF(16) *STATUS*J 

COMMON/*ADC1 / 1, ADC (16) 

V C0MM0N/*DAC1/  1 «DAC (16)  '*  , 

IERR=0  >;  . ; ’■ 

-'DO  1 1 = 1*16.  ; ,v 

)■  DAC(I)=OUT  > 1 / ■ 

1 BACK  ( I ) =ADC  ( I ) ‘ f ...  , v 

C COMPUTE  DIFFERENCE  IN  ADC  AND  DAC  VALUES  1 

C • V . • • • I 

• DO  2 1=1,16  . ■ 

DIFF  ( I ) =ABS (OUT) -ABS (BACK (1 ) ) 

IF(ABS(DIFF(I) ) .GT.BIGDIFF) IERR=1 

2 CONTINUE  • ' •-  . 

IF(IERR*EO.O)CALL  -SIMHOLD 

. STATUS=STAT(J)  :■  ■;  - ^ 

LOOP=LOOP*T  : 
f ' STATUS=STAT(J) 

* IF (LOOP*EO* 10) CALL  SIMHOLD 

CALL  SIMIOLE  ' 

- " END  ■ 1 . , 

0000000000000000000000  v 
RTREE  TRALGIO (0) *SUB1 ( 1 ) 

GLOBAL  'INTCOM 

END'  - f.W  ' i ' • v ' ' 

• oooooooooooooooooooooo 

000000000000000000000$  1 > . 
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